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SECTION  I 
INTRODUCTION 


The  motivation  for  this  report  was  to  develop  a  means  for 
editing  viscoelastic  material  data,  predicting  values  for 
parameters  of  equations  that  closely  describe  the  material  data, 
plotting  the  material  data,  and  storing  the  material  data  and 
parameters  by  use  of  a  fast,  dedicated  digital  computer  in  an 
interactive  mode  of  operation.  The  availability  of  a  computer 
program  to  accomplish  the  above  operations  would  greatly  decrease 
the  time  needed  to  characterize  the  viscoelastic  material  data. 

Up  to  this  point  in  time,  the  above-mentioned  material  data 
operations  had  been  accomplished  manually,  often  taking  many 
hours  of  calculations  and  plotting  to  represent  the  material  data. 

A  computer  program  could  and  would  eliminate  the  need  for  the 
manual  operations  of  calculating  and  plotting  the  data,  leaving 
only  the  manual  operations  of  providing  input  to  the  program  and 
operation  of  the  computer. 

The  analytical  representations  of  viscoelastic  material 
properties  have  evolved  over  many  years.  The  first  representation1 
displays  graphically  the  data  of  the  modulus  and  loss  factor 
properties  on  a  log-log  scale  with  an  independent  variable  called 
the  reduced  frequency.  This  variable  combines  the  temperature 
into  a  factor  that  multiplies  the  frequency,  the  factor  depending 
on  the  difference  between  the  temperature  and  a  suitably  chosen 
reference  temperature.  A  later  representation  displays  these 
properties  graphically  on  a  log-log  scale  with  an  independent 
variable  called  the  reduced  temperature.  The  reference  temperature 
is  unique  for  each  material  and  has  been  observed  experimentally 
to  be  approximately  equal  to  the  temperature  of  peak  damping 
for  many  materials. 

3 

Another  representation  displays  the  data  as  in  the  previous 
methods,  and  fits  the  properties  data  by  use  of  representative 
curves.  This  method  was  the  experimental  data  to  determine  a 
continuous  analytical  representation  for  the  material  properties. 
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This  report  will  include  descriptions  of  the  detailed  items 
that  are  features  of  the  computer  program.  Also  included  are  the 
equations  used  to  curve  fit  the  material  data  and  a  description  of 
the  curve-fitting  routine.  The  plotting  capabilities  will  be 
described,  a  user  manual  is  provided,  and  output  of  the  program 
is  documented  with  notes  to  indicate  the  meaning  of  the  input  and 
output  to  the  user.  The  listing  of  the  program  which  performs 
the  above  operations  on  the  viscoelastic  material  data  is  included 
as  an  appendix. 
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SECTION  II 


DEVELOPMENT  OF  THEORY 

In  order  to  facilitate  more  widespread  utilization  of 
viscoelastic  damping  material,  it  is  desirable  that  an  efficient, 
systematic,  and  standardized  procedure  be  established  to  process, 
store  ,  and  retrieve  the  material  properties.  The  approach 
adopted  is  to  use  equations  which  have  been  selected  to  best 
resemble  the  viscoelastic  material  properties.  A  minimization  of 
the  least  squares  will  yield  a  fit  to  the  material  data  for  the 
best  developed  equations  by  determining  the  values  of  parameters 
in  the  equations.  The  basis  for  data  processing,  storage,  and 
retrieval  is  the  use  of  these  parameters,  representing  the  best 
fit  of  the  material  data,  to  aid  in  efficient,  systematic,  and 
standardized  procedures  for  representing  the  viscoelastic  material 
properties  data. 

As  a  precursor  to  the  above  fit  of  the  equations,  an  editing 
procedure  has  been  developed  to  assist  in  the  sorting  of  data. 

The  material  data  can  be  expected  to  include  observations  which 
are  not  consistent  with  the  remaining  data,  and,  to  delete  this 
data  from  consideration,  an  editing  procedure  has  been  implemented. 
After  the  data  is  edited,  the  equations  are  fitted  to  the  remain¬ 
ing  material  data  for  further  processing. 

This  section  consists  of  subsections  representing  the  editing 
procedure,  the  empirical  shift  parameter  equation,  the  empirical 
modulus  equation,  and  the  empirical  loss  factor  equation. 

2.1  EDITING  PROCEDURE 

The  editing  procedure  consists  of  a  list  of  the  material  data 
and  several  plots.  The  plots  are  of  observed  frequency 
and  observed  loss  factor  versus  observed  temperature. 


The  loss  factor  versus  temperature  plot  is  important 
because  the  data  will  have  the  general  shape  of  an  inverted 
parabola,  and  a  smooth  curve  can  usually  be  fitted  through  the 
data.  From  this  smooth  curve,  the  peak  can  be  determined  and 
this  point  becomes  the  reference  temperature  (tQ)  used  to  "shift," 
and  minimize  the  scatter  of  the  modulus  and  loss  factor  data. 

The  frequency  versus  temperature  plot  is  important  because 

the  data  will  generally  lie  along  nearly  parallel  curves  for  each 
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resonant  mode  observed,  at  least  for  resonant  beam  type  * =sts  . 

The  plots,  along  with  a  list  of  the  material  data,  are  used 
to  edit  the  data.  If  ?.  data  point  does  not  lie  close  to  the 
general  parabolic  shape  of  the  loss  factor  versus  temperature 
curve  or  is  not  along  the  curves  of  frequency  versus  temperature, 
the  data  point  should  be  rejected  for  the  curve  fit,  since  the 
observation  generating  that  data  point  is  not  consistent  with  the 
remaining  observations. 

The  editing  procedure,  therefore,  is  to  eliminate  any  grossly 
inconsistent  data  points  or  observations  and  produce  a  reference 
temperature  (tQ)  for  the  shift  parameter  equation. 

2.2  EMPIRICAL  SHIFT  PARAMETER  EQUATION 

The  modulus  (M)  and  loss  factor  (n)  of  a  viscoelastic  material 
are  functions  of  two  independent  variables,  namely  temperature  (t) 
and  frequency  (f ) .  In  order  to  plot  the  modulus  and  loss  factor 
in  two  dimensions,  these  properties  must  be  functions  of  one 
variable,  where  this  one  variable  contains  the  two  independent 
variables  mentioned  above.  To  accomplish  this,  a  shift  parameter 
(aT) ,  which  is  a  function  of  the  temperature,  has  been  empirically 
developed1.  This  shift  parameter,  when  multiplied  by  the  frequency, 
gives  the  variable  of  reduced  frequency  (fr)  of  which  the  modulus 
and  loss  factor  are  functions. 

The  shift  parameter  equation  is  used  to  reduce  the  two 
dimensions  f  and  t  to  one  dimension  fat>  The  procedure  reduces 
the  scatter  of  the  dependent  modulus  and  loss  factor  equations 


I 

! 


which  have  been  reduced  from  surfaces  in  three  dimensions  with 
two  independent  variables,  to  curves  in  two  dimensions  with  one 

variable  (f  ) .  A  universal  shift  parameter  equation  has  been 
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empirically  determined  by  comparing  the  results  of  many  materials  . 

It  has  been  determined  that  the  shift  parameter  equation  is 
dependent  only  on  temperature  and  is  relative  to  a  reference  tem¬ 
perature  (tQ)  for  the  material,  which  reduces  the  scatter  of  the 
modulus  and  loss  factor  data  to  fit  a  smooth  curve  to  the  data. 

Based  on  rationale  presented  in  Reference  3,  the  "universal" 
shift  parameter  equation  to  be  used  for  the  processing  of  experi¬ 
mental  viscoelastic  material  data  is  as  follows: 


logiO (aT) 


-12(t-tQ) 


525 

.8 


in 


where  the  indicated  temperature  is  in  degrees  Centrigrade.  The 
reduced  frequency  is  computed  as  follows: 


therefore: 


fr  =  f-10 


-12 (t-t  )  1 

o 

rrl  +  (t“V 


or,  in  terms  of  the  log: 

12 (t-t  ) 

log10  (tr>  -  log10(f)  -  “S25  J 


rrs  * 


(2) 


This  equation  therefore  reduces  the  two  independent  variable  problem 
to  a  one  variable  problem  for  the  modulus  and  loss  factor  curves. 
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2.3  EMPIRICAL  MODULUS  EQUATION 


The  general  shape  of  the  log  modulus  versus  log  reduced 
frequency  curve  has  horizontal  asymptotes  about  which  it  is 
bounded.  The  curve  is  sigmoidal,  monotonic,  and  an  odd  function 
about  offsets.  The  equation  considered  in  this  report  was  derived 
in  Reference  3  and  is  of  the  general  form: 


log10(M) 


lo^10 (Mi} 
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1  + 


rom 
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(3) 


This  curve  is  asymptotic  to  M  for  low  values  of  reduced  frequency 

X / 

( f  ) ,  passes  through  the  offset  point  M  at  f  with  a  slope 
r  r  rom  rom 

governed  by  the  value  of  N,  and  has  an  asymptote  for  high  values 

of  reduced  frequency.  This  curve,  therefore,  is  an  odd  function 

about  the  point  (f  _ ,  M _ )  in  log-log  coordinates. 

rom  rom  3 

A  general  least  squares  fit  of  the  above  curve  to  the  material 

data  will  yield  parameter  values  for  f  ,  M  ,  N,  and  M.  that 

rom  rom  It 

will  best  describe  the  material  data,  where  M  is  the  observed 

modulus  and  f  is  the  observed  reduced  frequency  presented  in 

r  2 
Section  2.2.  The  indicated  modulus  values  are  in  N/m  ,  the  indicated 

reduced  frequency  values  are  as  presented  in  Section  2.2,  and  the 

slope  parameter  value  at  the  offset  is  dimensionless  in  the  presented 

modulus  curve. 


This  curve  can  be  used  to  represent  any  type  of  modulus.  In 
the  case  of  shear  modulus,  simply  replace  the  M's  in  Equation  3 
with  G's,  indicating  shear  modulus.  In  the  case  of  Young's  modulus, 
simply  replace  the  M's  in  Equation  3  with  E's  indicating  Young's 
modulus .  The  equation  of  the  modulus  curve  presented  can  be  used 
to  fit  the  material  data  to  a  smooth  curve,  and  can  also  be 
represented  on  a  graph  of  the  viscoelastic  material  data.  The 
parameter  values  determined  by  the  least  squares  fit  can 


be  stored  on  the  computer,  and  the  equation,  along  with  the 
parameters,  can  be  used  at  a  future  time  to  predict  the  modulus  of 
the  material  for  which  the  parameters  are  defined.  Therefore, 
the  storage  and  retrieval  of  data  from  the  computer  is  important 
for  further  studies  of  viscoelastic  material  properties. 

2.4  EMPIRICAL  LOSS  FACTOR  EQUATION 

The  general  shape  of  the  log  loss  factor  versus  log  reduced 
frequency  curve  is  hyperbolic  in  nature.  The  general  curve  is 
an  offset  symmetric  hyperbola  with  straight  line  asymptotes 
added  to  the  equation. 

The  equation  considered  was  derived  in  Reference  3  and  is 
of  the  general  form: 


log10(n>  -  log10(nfrot)  ♦  §  <Vsh>l  lo*l0 
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This  curve  has  a  peak  of  nfro4  at  f  A  with  a  peak  shape  defined 
by  C,  is  asymptotic  to  a  line  of  slope  S^  for  low  values  of  reduced 
frequency  (fr) ,  and  is  asymptotic  to  a  line  of  slope  for  high 
values  of  reduced  frequency. 


A  general  least  squares  fit  of  the  above  curve  to  the 
material  data  will  yield  parameter  values  for  gj,'  Sh' 


f  A,  and  C  that  will  best  describe  the  material  data,  where  n 


is  the  observed  loss  factor  and  fr  is  the  observed  reduced  frequency 


present  in  Section  2.2.  The  indicated  reduced  frequency  values 
are  as  presented  in  Section  2.2,  and  the  remaining  values  are 
dimensionless  in  the  presented  loss  factor  curve. 

The  equation  of  the  loss  factor  curve  can  be  used  to  fit 
the  material  data  to  a  smooth  curve,  and  can  also  be  represented 
on  a  graph  of  the  viscoelastic  material  data.  The 


r, 


parameter  values  determined  by  the  least  squares  fit  can  be  stored 
on  the  computer,  and  the  equation  along  with  the  parameters  can  be 
used  at  a  future  time  to  predict  the  loss  factor  of  the  material 
for  which  the  parameters  are  defined. 

2 . 5  PARAMETER  CONSTRAINTS 
3 

A  study  was  performed  to  examine  the  effect  of  changing 
parameters  and  applying  the  equations  to  a  limited  span  of  obser¬ 
vations,  and  it  was  found  that  constraints  had  to  be  imposed  on 
the  parameters  for  the  modulus  and  loss  factor  curves.  The  study 
concluded  that  the  least  squares  fit  to  the  material  data  gave  an 
excellent  fit  to  the  span  data,  but  that  the  parameter  values 
obtained  were  not  realistic  for  the  physical  interpretation  of  the 
equations.  Also,  the  study  concluded  that  parameters  should  be 
held  constant  is  possible  and  also  may  be  constrained  for  the 
physical  interpretation  of  the  curves. 
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Based  on  the  least  squares  fit  in  the  study  ,  the  follow¬ 
ing  constraints  for  the  given  parameters  are  shown  in  general: 
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(5) 


By  imposing  the  above  constraints  and  allowing  parameters  to  be 
held  constant,  the  curve  fits  can  now  take  on  physical  meaning. 


To  obtain  the  "best"  possible  fit  to  the  material  data, 
a  wide  range  of  observations  must  be  taken  to  ensure  the 
material  data  can  be  characterized  by  the  empirical  curves 
presented.  The  studyJ  shows  that  the  viscoelastic  material 
data  can  be  characterized  by  the  equations,  and  with  some  degree 
of  certainty  these  equations,  along  with  the  desired  material 
parameters,  can  be  used  to  predict  the  viscoelastic  material 
properties. 
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SECTION  III 
PROGRAM  DEVELOPMENT 


The  computer  program  wac  developed  as  a  result  of  discussions 
with  Air  Force  Materials  Laboratory  (AFML)  and  Air  Force  Flight 
Dynamics  Laboratory  (AFFDL)  personnel  concerned  with  the  charac¬ 
terization  of  material  properties  and  usage  of  viscoelastic 
materials  in  damping  treatments.  The  type  of  program  developed 
was  beneficial  for  the  physical  representation  of  the  viscoelastic 
material  data  and  utilizations  by  the  intended  users.  The 
development  included  the  general  flow  for  the  program,  the  use  of 
available  equipment,  and  the  use  of  a  curve  fit  routine  to  fit 
the  data  to  the  nonlinear  curves  presented  in  Section  II.  These 
items  were  the  basis  for  the  program  development. 

3.1  GENERAL  PROGRAM  FLOW 

A  general  flow  of  the  computer  program  was  designed,  as 
presented  in  Figure  1,  with  the  ability  to  read  in  viscoelastic 
material  property  data,  edit  the  data  to  eliminate  any  inconsistent 
data,  obtain  a  reference  temperature  that  reduces  the  scatter  of 
the  data  about  the  given  curves,  fit  the  modulus  data  to  a  charac¬ 
teristic  curve,  and  finally,  plot  the  edited  viscoelastic  material 
property  data  on  the  desired  grid  with  corresponding  labels. 

As  the  program  developed  and  after  further  discussions  with 
AFML  and  AFFDL  personnel,  extensions  to  the  basic  program  were 
instituted  to  make  the  program  more  desirable  for  the  user.  An 
option  was  made  available  to  store  the  data  and  parameters  of  the 
curves  within  the  computer  facility. 

Since  a  store  option  was  included,  a  retrieve  option  was  added 
so  the  user  can  store  the  data  at  one  time  and  execute  the  program 
at  a  future  time,  retrieve  the  data,  and  continue  the  analysis  of 
the  data.  Also  within  the  edit  routine,  an  added  option  was 
included  so  the  user  could  retrieve  stored  data,  add  more  data  for 
the  analysis,  and  have  the  option  to  store  the  expanded  data  set. 
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Also,  at  this  point  a  decision  was  made  to  perform  all  calculations 
and  final  plots  in  metric  units,  therefore,  an  option  was  included 
to  allow  the  user  to  read  in  new  data  in  either  metric  or  English 
units,  and  to  add  new  data  in  the  edit  routine,  in  metric  or  English 
units.  A  more  detailed  general  flow  of  the  program  is  shown  in 
Figure  2. 

Further  discussions  with  AFML  and  AFFDL  personnel  have  led 
to  more  additions  to  the  basic  program.  These  additions  do  not 
affect  the  flow  of  the  program  presented.  These  additions  are 
basically  items  that  interface  with  the  user-program  interaction 
within  each  block  in  the  figure. 

3.2  AVAILABLE  EQUIPMENT 

The  program  development  was  limited  to  the  equipment 
available  for  implementation  and  usage  of  the  computer  program. 

The  computer  program  was  coded  in  FORTRAN  IV  for  the  HP-2100 
minicomputer  under  the  DOS-Ill  operating  system.  This  minicomputer 
is  available  at  the  AFML  where  the  implementation  and  usage  of 
the  program  is  to  be  available  for  AFML  personnel.  The  DOS-Ill 
operating  system,  with  a  FORTRAN  IV  compiler,  allows  the  storage 
of  data  on  files  resident  on  a  disk.  This  feature  is  useful  in 
the  program  and  in  the  program  development  by  implementing  disk 
storage  of  data  by  the  program.  The  disk  storage  capability  also 
allows  the  user  to  read  and  write  data  to  files,  and  therefore 
to  build  a  thorough  data  base  of  viscoelastic  material  property 
data. 

Also  available  at  the  AFML  is  a  Tektronix  graphics  screen 
terminal,  which  is  needed  to  accomplish  the  required  computer 
generated  graphics  of  the  program.  This  terminal  is  connected 
directly  to  the  HP-2100  computer  mentioned  above.  The  terminal 
has  a  hard-copy  unit  which  produces  a  copy,  on  a  sheet  of  paper, 
of  all  information  projected  onto  the  screen  of  the  terminal. 

This  equipment  is  valuable  for  generating  the  plots  of  the  material 
data,  and  has  the  ability  to  obtain  paper  copies  of  the  screen  for 
reference  at  a  future  time. 
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Figure  2.  Detailed  Program  Diagram. 


The  software  required  for  the  implementation  of  this  program 
is  the  DOS-Ill  operating  system  which  allows  storage  of  data, 
the  FORTRAN  IV  compiler,  and  a  package  of  utility  routines  used 
to  project  the  graphics  or  plots  onto  the  screen  of  the  terminal. 
These  utility  routines  are  commonly  known  as  the  PLOT-10  software 
package  for  the  Tektronix  hardware®  These  routines  allow  the 
user  to  display  any  information  on  any  part  of  the  screen.  This 
is  a  very  powerful  package  for  projecting  graphics  and  alpha- 
numerics  on  the  screen  of  the  terminal. 

The  decision  was  made  to  use  the  above  mentioned  equipment 
and  software.  Since  AFML  has  the  facility  indicated,  the 
program  development  and  implementation  was  completed  at  AFML, 
for  the  purpose  of  program  use  by  AFML  personnel  responsible 
for  the  characterization  of  viscoelastic  material  property  data. 

3.3  TYPE  OF  MINIMIZATION 

The  program  development  also  led  to  the  examination  of  curve 
fitting  routines  for  the  modulus  and  loss  factor  curves  presented 
in  Section  II.  The  decision  was  made  to  use  the  minimization 
technique  known  as  least-squares  to  determine  the  parameters  that 
most  closely  describe  these  nonlinear  curves  of  the  data.  In  this 
case,  the  least-squares  technique  leads  to  the  minimization  of  a 
nonlinear  function,  being  the  sum  of  the  square  deviations  of  the 
data  about  the  defining  curve.  A  direct  search  approach  was  used 
to  scan  the  parameter  space,  within  which  the  nonlinear  curve  is 
defined,  and  determine  a  point  in  the  parameter  space  that  minimizes 
the  nonlinear  function.  The  direct  search  approach,  as  applied  to 
the  computer  program,  implements  some  type  of  prescribed  search 
pattern  to  arrive  at  the  minimum  of  the  nonlinear  function  being 
considered.  Once  this  point  has  been  found,  the  prescribed  search 
pattern  determines  that  any  movement  in  the  parameter  space  by  the 
direct  search  approach  will  only  yield  an  increase  in  the  value  of 
the  nonlinear  function  being  minimized.  Therefore,  this  point 
must  be  the  minimum  of  the  nonlinear  function  within  a  region  of 
the  minimizing  point  in  the  parameter  space. 


This  type  of  minimization  procedure  was  chosen  because  it 
does  search  the  parameter  space  directly  to  find  the  minimum  value 
of  the  function  to  be  minimized.  Also,  this  type  required  a  much 
smaller  amount  of  coding  than  other  techniques,  thus  allowing 
more  coding  for  other  parts  of  the  overall  program.  Also,  this 
type  is  relatively  fast  for  small  dimensionality  of  the  parameter 
space  being  considered. 
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SECTION  IV 
PROGRAM  DESCRIPTION 

This  section  contains  an  overall  general  description  of  the 
program  implemented  to  represent  the  viscoelastic  material  proper¬ 
ties  data.  Several  items  that  will  be  discussed  are  the  interactive 
nature  of  the  program  including  inputs  and  outputs,  a  description 
of  the  minimization  technique,  a  description  of  the  graphics 
package,  a  description  of  the  structure  of  the  data  files  for  data 
storage,  and  a  description  of  the  initialization  of  the  parameter 
values  for  the  modulus  and  loss  factor  curves  presented  in 
Section  2. 

4 . 1  INTERACTIVE  FEATURES 

The  computer  program  was  developed  to  accommodate  the  inter¬ 
active  features  requested  by  AFML  and  AFFDL  personnel.  The  user 
manual,  presented  in  Section  6,  presents  all  the  possible  input 
and  output  features  available  in  the  program  which  require  user 
interaction.  The  interactive  features  provided  are  of  three 
basic  types. 

The  decision  type  involves  alphanumeric  input  by  the  user  for 
decision  type  of  responses  for  the  program.  These  input  types  are 
preceded  by  the  output  of  a  question,  for  which  the  user  inputs 
the  correct  alphanumeric  character  for  the  response  to  the  question. 
The  program  then  interprets  this  decision  and  proceeds  to  the  next 
sequential  segment  of  the  program  based  on  the  response  decision. 
These  decision  alternatives  are  presented  in  Section  6,  the  User's 
Manual,  for  each  question  output  by  the  program  which  requires  a 
decision  type  of  response  by  the  user. 

The  numeric  type  of  input  involves  numeric  values  to  be 
entered  by  the  user.  These  inputs  are  preceded  by  an  output,  for 
which  the  user  inputs  the  correct  numeric  values  for  the  response 
to  the  output.  This  output,  in  general,  implies  that  the  user  is 
to  enter  data  or  values  for  the  following  input.  This  input  must 
be  in  the  form  of  a  number,  thus  indicating  a  numeric  value  being 
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entered  to  the  program.  These  numeric-type  inputs  are  presented 
in  Section  6,  the  User's  Manual,  for  each  output  by  the  program 
which  requires  a  numeric  type  of  input  or  response  by  the  user. 

The  control  type  input  involves  the  entering  of  control 
instructions  by  the  user.  These  control  instructions  are  instruc¬ 
tions  to  the  operating  system  for  the  HP-2100  minicomputer.  The 
user  inputs  a  control  instruction  to  initiate  execution  of  the 
program.  After  execution  begins,  the  program  will  periodically 
suspend  and  pause  execution  for  the  user  to  obtain  copies  of  the 
terminal  screen.  The  user  then  inputs  a  control  instruction  to 
continue  execution  of  the  program  or  discontinue  execution  of 
the  program,  causing  an  abnormal  termination  of  the  program. 

These  control  inputs  are  presented  in  Section  6,  User's  Manual. 

4.2  DATA  FILE  STRUCTURE 

The  computer  program  development  led  to  the  discussion  of 
data  storage  with  AFML  and  AFFDL  personnel.  The  discussion 
led  to  three  data  files  for  data  storage  by  the  program.  A  general 
description  of  these  files  and  type  of  data  acquisition  is 
discussed. 

The  description  of  these  data  files  is  rather  simple.  The 
program  development  led  to  the  use  of  three  files  to  store  data. 

The  first  file  is  a  general  file  which  stores  the  material  name, 
data  sources,  and  other  items  (i.e.,  number  of  data  points,  material 
code  number,  data  type,  etc.)  for  the  material  on  the  designated 
material  code  number.  This  number  represents  the  storage  location 
in  the  data  file  for  the  above-mentioned  data.  This  file  is  used 
to  indicate  if  any  material  data  is  stored  on  the  designated 
material  code  (see  Section  6) .  The  second  file  is  the  data  file 
to  store  the  reduced  and  edited  material  data.  This  data  file 
can  be  accessed  through  the  first  file  and  the  material  data  can 
be  found  relative  to  the  material  code  number.  This,  again, 
represents  the  storage  location  in  the  data  file  for  the  reduced 
and  edited  material  data. 
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The  third  file  is  the  data  file  to  store  the  parameter  values 
determined  by  the  curve  fits  of  this  program.  This  data  file  can 
be  accessed  through  the  first  file  and  the  parameter  data  can  be 
found  relative  to  the  material  code  number.  This  again  represents 
the  storage  location  in  the  data  file  for  the  parameter  values. 

This  third  file  cannot  be  retrieved  directly  by  the  user  of  this 
program  as  the  other  two  can  (Section  6) .  This  file  is  the  product 
of  the  computer  program  and  will  be  used  at  a  future  time  in 
another  program,  which  will  use  the  parameter  values  for  a  material 
to  predict  the  material  properties.  As  can  be  seen,  these  three  files 
are  combined  by  use  of  the  material  code  number  of  the  first  file 
to  completely  store  the  data  for  this  program  and  subsequent 
follow-on  programs.  This  type  of  data  file  structure  eliminates 
the  need  for  the  user  to  enter  the  entire  data  at  the  terminal 
at  each  execution  of  the  program.  This  structure  also  allows  the 
building  of  data  sets  for  each  material  in  the  edit  sequence 
of  the  program,  and  then  store  the  expanded  data  set  for  the 
material.  This  structure  also  allows  re-editing  and  re-fitting 
of  the  data  to  the  curves  until  a  representative  graph  of  the 
material  data  has  been  achieved,  without  the  need  of  re-entering 
the  data  each  time.  Finally,  this  structure  allows  the  implementa¬ 
tion  of  other  programs  to  store  data  on  these  files  for  use  by 
this  program,  and  to  retrieve  data  from  these  files  for  use  by 
subsequent  analysis-type  programs. 

The  above  mentioned  data  files  are  resident  on  the  HP-2100 
minicomputer  disk  at  the  AFML,  along  with  the  computer  program. 

The  operating  system  allows  the  programmer  to  implement,  in  the 
program,  direct  storage  and  retrieval  of  data  on  the  disk.  This 
data  must  have  specified  file  names  and  be  able  to  locate  storage 
areas,  in  the  files,  relative  to  the  beginning  of  the  files.  For 
this  computer  program,  the  three  files  are  implemented.  By  use 
of  a  utility  routine  of  the  operating  system  callable  by  FORTRAN  IV, 
the  program  can  store  or  retrieve  (write  or  read)  data  from  the 
files  mentioned  and  from  the  storage  location  indicated  by  the 
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material  code  for  each  file.  This  capability  of  the  operating 
system  assisted  in  the  development  of  this  computer  program  and 
aids  in  the  storage  and  retrieval  of  material  data.  Thus, 
having  stored  a  massive  amount  of  material  properties  data,  these 
files  can  be  used  for  curve  fitting,  further  analysis,  and 
predictions  of  the  viscoelastic  material  properties. 

4 . 3  GRAPHICS  PACKAGE 

The  hardware  to  accomplish  the  required  graphics  is  the 
TEKTRONIX  4014-1  graphics  screen  terminal.  This  terminal  is 
well  adapted  for  graphics  because  it  has  what  is  commonly  known  as 
a  vector  mode  for  the  terminal.  In  the  vector  mode  the  programmer 
has  the  ability  to  address  individual  points  on  the  screen,  thus 
having  the  ability  to  draw  lines  and  have  complete  control  over 
the  addressibility  of  the  entire  screen.  In  the  alphanumeric 
mode,  the  programmer  can  use  the  FORTRAN  IV  input/output  to  display, 
on  the  screen,  the  input  and  output  required  by  the  program.  This 
hardware  is  connected  directly  to  the  HP-2100  mini-computer  on 
which  this  program  is  to  be  implemented. 

To  accomplish  graphics  on  the  terminal,  the  TEKTRONIX  Corp. 
has  supplied  a  package  of  subprograms  which  can  be  used  by  the 
programmer  to  generate  the  graphics  and  thus  project  the  graphics 
on  the  screen  of  the  terminal  for  display  to  the  user.  The 
package  supplied  involves  complete  terminal  control  subprograms 
which  control  the  status  of  the  terminal  (i.e.,  vector  mode  or 
alphanumeric  mode)  so  the  programmer  has  complete  control  of  the 
terminal  at  all  times  in  the  program.  This  type  of  programmer 
control  allows  the  coexistence  of  graphics  and  FORTRAN  input/output 
on  the  screen  for  display  to  the  user.  The  subprograms  implemented 
for  this  program  are  not  described  in  the  subprogram  description 
(Section  7) .  These  subprograms  are  described  in  Reference  6,  which 
includes  a  description  of  each  subprogram,  what  the  subprogram 
accomplishes  for  the  terminal,  and  how  the  subprogram  controls  the 
display  on  the  screen  of  the  terminal.  This  graphics  package. 
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along  with  the  hardware,  allows  the  programmer  to  display  any 
required  graphics  and/or  input/output  required  for  the  implemen¬ 
tation  of  the  computer  program  to  generate  a  representation  of 
the  viscoelastic  material  properties. 

4.4  MINIMIZATION  TECHNIQUE 

The  nonlinear  function  described  for  the  minimization  is 
the  sum  of  square  deviations  or,  as  commonly  known,  least  squares 
fit  of  data  to  a  given  curve.  The  type  of  minimization  imple¬ 
mented  in  the  computer  program  is  the  direct  search  method.  In 
this  method,  the  algorithm  searches  the  parameter  space  directly 
to  determine  the  minimum  of  the  function.  To  find  the  values  of 
the  parameters  in  the  parameter  space,  this  method  uses  some 
predetermined  search  pattern  to  locate  the  minimum  of  the 
function  in  the  given  parameter  space.  The  direct  search  technique 
uses  what  is  commonly  known  as  a  SIMPLEX  in  the  search  pattern. 

A  SIMPLEX  is  a  collection  of  points  in  the  parameter  space, 
for  which  each  point  is  called  a  vertex  and  is  represented  by 
the  function  value  and  parameter  values.  The  direct  search  method 
replaces  each  vertex  with  another  vertex  which  yields  a  smaller 
functional  value.  This  scheme  continues  until  all  the  vertices 
converge  to  one  point  which  should  be  the  minimum  functional  value. 
This  point,  therefore,  represents  the  function  and  parameter  values 
that  yield  a  minimum  for  the  function  within  a  given  neighborhood 
of  the  point  in  the  parameter  space.  Any  search  efforts  beyond 
this  point  will  only  yield  a  larger  value  for  the  function  in  the 
parameter  space.  These  parameter  values  of  the  point,  therefore, 
represent  the  parameters  that  force  the  function  to  be  a  minimum. 

In  this  program  the  function  is  again  the  least  squares  fit  of  data 
to  desired  curves  (Section  2) ,  and  the  parameter  values  are  those 
values  for  the  parameters  of  the  curve  that  cause  the  function  to 
be  a  minimum.  This  direct  search  technique,  therefore,  has  found 
parameter  values  for  the  desired  curves  that  most  closely  represent 
the  viscoelastic  material  properties. 
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4.5  INITIALIZATION  OF  PARAMETER  VALUES 


The  minimization  technique  requires  initial  parameter 
values  for  the  direct  search  procedure.  These  initial  values 
are  used  to  provide  a  starting  point  for  the  search  and  are  also 
used  to  describe  the  search  pattern  required  for  the  direct 
search  method. 

For  the  modulus  curve,  the  initial  parameter  values  are 

assigned  in  the  following  way.  The  f  parameter  value  is  found 

by  taking  the  square  root  of  the  product  of  the  minimum  and 

maximum  observed  reduced  frequency  values  from  the  data.  The 

Mrom  Parameter  initial  value  is  found  by  fitting  a  second-degree 

polynomial  to  the  log  modulus  versus  log  reduced  frequency  data, 

substituting  the  value  of  the  log  f  above  to  determine  the 

rom 

value  of  log  Mrom#  and  then  looking  up  the  anti-log  to  find  the 
value  for  Mrom.  The  M  parameter  initial  value  is  found  by 
taking  one-half  of  the  minimum  observed  modulus  value  from  the 
data.  The  N  parameter  initial  value  is  assigned  to  .4. 

For  the  loss  factor  curve,  the  initial  parameter  values 
are  assigned  in  the  following  way.  The  f  .  parameter  initial 
value  is  found  by  fitting  a  second-degree  polynomial  to  the  log 
loss  factor  versus  log  reduced  frequency  data,  determining  the 
location  of  the  peak  for  the  polynomial,  and  taking  the  anti-log 
of  the  log  reduced  frequency  value  to  assign  the  value  of  fro£- 
The  1jr0£  parameter  initial  value  is  found  by  determining  the 
value  of  the  peak  for  the  polynomial  fit,  and  taking  the  anti-log 
of  the  log  loss  factor  to  assign  the  value  of  n^rQ^.  The 
parameter  initial  value  is  assigned  to  .35.  The  parameter 
initial  value  is  assigned  to  -.45.  The  C  parameter  initial  value 
is  assigned  to  4.0. 

The  initial  parameter  values  are  assigned  by  the  above 
methods.  The  search  technique  uses  these  initial  values  to 
determine  the  search  procedure  to  implement  in  the  parameter  space. 
The  search  technique  then  proceeds,  from  the  initial  values,  to 
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search  the  parameter  space  and  determine  the  parameter  values 
that  will  most  closely  describe  the  desired  curves  for  the 
viscoelastic  material  properties  data. 


SECTION  V 

SUBPROGRAM  DESCRIPTION 


This  section  includes  a  general  description  of  the  routines 
which  are  combined  together  to  form  the  computer  program  which 
generates  the  representation  of  the  viscoelastic  material  proper¬ 
ties  data.  Many  of  the  subprogram  routines  provide  user  inter¬ 
action  and  input  to  control  the  representation  of  the  viscoelastic 
material  properties.  This  section  does  not  include  a  description 
of  the  graphic  routines  needed  for  the  overall  computer  program. 

A  description  of  these  routines  can  be  found  in  Reference  6 . 

5.1  PROGRAM  VIBR2 

This  routine  is  the  main  segment  of  the  computer  program. 

When  execution  of  the  program  begins,  this  segment  initializes 
the  graphics  package  for  further  use  by  the  program.  This  segment 
makes  several  calls  to  routines  in  the  graphics  package  to  set  the 
appropriate  values  to  be  used  by  the  graphics  package.  This  seg¬ 
ment  first  calls  the  read-in  and  retrieval  routine  to  make  avail¬ 
able  the  data  to  be  represented.  This  segment  then  inquires  about 
editing  the  listed  data.  If  so,  this  segment  calls  the  editing 
routine.  In  either  case,  this  segment  then  calls  the  routine  to 
fit  the  loss  factor  data  to  a  second  degree  polynomial  curve  of 
the  temperature  data.  This  segment  then  calls  the  graphics 
routine  which  generates  two  graphs  and  lists  the  data.  This  seg¬ 
ment  again  inquires  about  editing  the  data.  If  so,  the  editing 
routine  is  called,  and  this  segment  proceeds  back  to  call  the 
routine  to  fit  the  loss  factor  versus  temperature  data.  If  not, 
this  segment  inquires  about  the  changing  of  the  reference  tempera¬ 
ture.  If  so,  this  segment  has  the  user  input  the  reference  tem¬ 
perature.  In  either  case  this  segment  calls  the  routine  to  fit  the 
modulus  data  to  the  specified  curve  presented  in  Section  2.  This 
segment  then  calls  the  routine  to  fit  the  loss  factor  data  to  the 
specified  curve  presented  in  Section  2.  This  segment  inquires  about 
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the  user  input  for  the  interval  of  temperature  lines  for  the  repre¬ 
sentative  graph.  This  segment  then  calls  the  graphing  routine 
which  generates  the  appropriate  grid,  graphs  the  data,  and  graphs 
the  curves  determined  by  the  above  curve  fits.  This  segment  then 
calls  the  routine  which  outputs  the  curves  used  in  the  curve  fits 
and  the  parameters  that  describe  the  curves  for  the  listed  data. 

This  segment  inquires  about  changing  the  reference  temperature  and/or 
parameters.  If  so,  this  segment  has  the  user  input  the  reference 
temperature  and  proceeds  back  to  call  the  modulus  curve  fitting 
routine.  If  not,  this  segment  inquires  about  the  storing  of  the 
data.  If  so,  this  segment  calls  the  storing  routine.  In  either 
case,  this  segment  inquires  about  the  use  of  the  program  to 
represent  another  set  of  data.  If  so,  this  segment  proceeds  back 
to  the  read-in  and  retrieval  routine  after  clearing  the  screen  of 
the  terminal  to  prepare  another  set  of  data  to  be  represented.  If 
not,  this  segment  terminates  the  graphics  package  and  terminates 
execution  of  the  program.  Control  is  then  returned  to  the 
operating  system. 

5.2  SUBROUTINE  REDIN 

This  segment  of  the  program  is  used  to  make  available  the 
data  to  be  represented.  This  segment  inquires  about  the  listing 
of  labels  for  previously  stored  data.  If  so,  this  segment,  by  use 
of  a  utility  routine  of  the  operating  system,  displays  the  labels 
of  the  desired  stored  data.  In  either  case,  this  segment  inquires 
about  the  reading  in  of  new  data.  If  so,  this  segment  asks  many 
questions  and  requires  much  input  from  the  user  to  characterize 
and  label  the  data.  The  user  has  the  option  to  enter  the  data  in 
metric  units  or  English  units.  The  user  then  inputs  the  data  in 
either  unit  and,  after  it  is  entered,  the  data  is  converted  to  metric 
units  because  all  calculations  and  computations  are  accomplished 
in  metric  units.  If  not,  this  segment  inquires  a  retrieval  code 
which  represents  data  that  has  already  been  stored  by  the  program 
or  by  other  programs.  This  segment  uses  a  utility  routine  of  the 
operating  system  to  retrieve  the  stored  data  from  the  files  of  the 
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computer  system.  In  either  case ,  this  segment  erases  the  screen 
and  redisplays  this  newly  read-in  data  or  retrieved  data  now 
output  in  metric  units.  This  segment  then  pauses  in  its  execution 
of  the  program,  so  that  the  user  at  his  option  can  obtain  copies 
of  the  terminal  screen  for  further  reference,  and  then  returns 
control  of  the  program  to  its  calling  routine. 

5.3  SUBROUTINE  FITTO 

This  segment  is  used  to  fit  a  second  degree  polynomial 
to  the  loss  factor  data  as  a  function  of  the  temperature,  by  use 
of  a  routine  to  fit  the  data.  This  segment  also  determines  the 
peak  of  the  polynomial  to  give  the  reference  temperature  of  the 
material  to  the  nearest  multiple  of  10 °C.  This  segment  then 
returns  control  to  the  calling  routine. 

5.4  SUBROUTINE  POLYF 

This  segment  is  used  to  fit  a  second  degree  polynomial 
to  any  data.  This  segment  is  included  because  the  program  needs, 
in  several  instances,  a  polynomial  fit  to  the  data.  This  segment 
then  returns  control  to  its  calling  routine. 

5.5  SUBROUTINE  PLTTO 

This  segment  of  the  program  is  used  to  set  the  specified 
grid  size  for  the  frequency  versus  temperature  and  loss  factor 
versus  temperature  graphs.  This  segment  clears  the  screen  of 
the  terminal,  sets  the  specified  grid  sizes,  and  then  calls  a 
plotting  routine  twice  to  generate  the  two  graphs,  mentioned  above, 
on  the  screen.  This  segment  then  prepares  the  terminal  for  further 
output  by  several  calls  to  routines  of  the  graphics  package.  This 
segment  then  outputs  the  data  that  is  represented  by  the  graphs 
displayed  on  the  screen,  and  returns  control  of  the  program  to  its 
calling  routine. 
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5.6 


SUBROUTINE  PLTT1 


This  segment  of  the  program  is  used  to  generate  a  graph. 
This  segment  generates  either  the  frequency  versus  temperature 
plot  or  the  loss  factor  versus  temperature  plot.  For  the  former, 
this  segment  outputs,  through  calls  to  the  routines  of  the  graphics 
package,  the  grid  of  the  plot,  the  labels  of  the  axes  for  the  grid, 
and  the  plotting  of  the  points,  along  with  an  identifying  number, 
of  the  frequency  versus  temperature  data.  For  the  latter,  this 
segment  outputs,  through  calls  to  routines  of  the  graphics  package, 
the  grid  of  the  plot,  the  labels  of  the  axes  for  the  grid,  the 
plotting  of  the  points,  along  with  an  identifying  number,  of  the 
loss  factor  versus  temperature  data,  a  curve  through  the  loss 
factor  versus  temperature  data  mentioned  above,  and  an  output  of 
the  reference  temperature  for  the  material  data.  This  segment 
then  returns  control  to  its  calling  routine. 

5 . 7  SUBROUTINE  EDITO 

This  segment  is  used  to  edit  the  material  data  which 
has  been  read  in  or  retrieved.  The  options  available  to  the 
user  are  the  ability  to  add  new  data,  delete  data,  or  end  the 
edit  sequence.  If  the  user  wants  to  add  new  data,  this  segment 
allows  adding  of  data  in  metric  or  English  units.  The  user  then 
inputs  the  new  data,  and  the  new  data  is  then  converted  to 
metric  units  after  input.  If  the  user  wants  to  delete  data,  this 
segment  requires  input,  by  the  user,  for  the  lines  or  point  numbers 
from  the  data  list  to  be  deleted,  These  points  are  deleted  and 
the  data  is  resequenced  to  indicate  that  the  data  has  been  entirely 
deleted  from  the  data  being  considered.  If  the  user  wants  to  end 
the  edit  sequence,  this  segment  terminates  and  returns  to  its 
calling  routine, 

5 . 8  SUBROUTINE  STORE 

This  segment  of  the  program  is  used  to  store  the  material 
label,  material  data,  and  parameters  of  the  curve  fit  (Section  2) 
on  specified  files.  These  files  are  set  to  be  data  files  for  this 
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program  and,  through  the  use  of  a  utility  routine  of  the  operating 
system,  this  data  can  be  stored  or  written  to  these  files  on  the 
disk  of  the  computer.  This  data  is  stored  for  future  reference 
and  possibly  future  analysis  needed  for  the  given  material.  This 
segment  outputs  a  message  that  the  data  has  been  stored  and  pauses 
execution  of  the  program  so  the  user  can  obtain  copies  of  the 
terminal  screen  for  further  reference,  and  returns  control  of  the 
program  to  its  calling  routine. 

5 . 9  SUBROUTINE  MODUS 

This  segment  of  the  program  is  used  to  fit  the  modulus  data 
to  the  curve  presented  in  Section  2.  This  segment  clears  the  screen 
of  the  terminal  and  outputs  a  message  to  the  user.  This  segment 
then  initializes  the  parameters  for  the  modulus  curve  fit  from  the 
given  data  with  the  use  of  the  polynomial  fit  routine.  This  seg¬ 
ment  then  inquires  about  the  constraints,  for  which  the  user 
enters  the  constraint  values  as  presented  in  Section  2.  This  seg¬ 
ment  then  calls  the  minimization  routine  to  determine  the  values  of 
the  parameters  for  the  modulus  curve  fit.  This  segment  calls  the 
routine  to  output  the  parameter  values  and  modulus  curve  fit  to  the 
terminal.  This  segment  then  inquires  if  the  user  wants  to  accept 
these  parameter  values.  If  he  does,  this  segment  returns  control 
to  its  calling  routine.  If  not,  this  segment  inquires  about  the 
initial  values  the  user  wants  to  input  for  the  parameters,  and 
the  ability  of  the  user  to  hold  parameters  constant  in  the  mini¬ 
mization  process.  This  segment  then  proceeds  back  to  call  the 
minimization  routine. 

5.10  SUBROUTINE  ETA 

This  segment  of  the  program  is  used  to  fit  the  loss  factor 
data  to  the  curve  presented  in  Section  2.  This  segment  clears  the 
screen  of  the  terminal  and  outputs  a  message  to  the  user.  This 
segment  then  initializes  the  parameters  for  the  loss  factor  curve 
fit  from  the  given  data  with  the  use  of  the  polynomial  fit  routine. 
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This  segment  then  inquires  about  the  constraints  for  which  the 
user  enters  the  constraint  values  as  presented  in  Section  2. 

This  segment  then  calls  the  minimization  routine  to  determine  the 
values  of  the  parameters  for  the  loss  factor  curve  fit.  This  seg¬ 
ment  calls  the  routine  to  output  the  parameter  values  and  loss 
factor  curve  fit  to  the  terminal.  This  segment  then  inquires  if 
the  user  wants  to  accept  these  parameter  values.  If  he  does,  this 
segment  returns  control  to  its  calling  routine.  If  not,  this  seg¬ 
ment  inquires  about  the  initial  values  the  user  wants  to  input  for 
the  parameters  and  the  ability  of  the  user  to  hold  parameters  con¬ 
stant  in  the  minimization  process.  This  segment  then  proceeds  back 
to  call  the  minimization  routine. 

5.11  SUBROUTINE  FCN 

This  segment  is  used  by  the  minimization  routine  to  obtain 
the  value  of  the  function  to  be  minimized,  namely,  the  sum  of 
square  deviations.  This  segment  uses  the  reference  temperature 
routine  to  obtain  the  value  of  the  independent  variable,  namely, 
reduced  frequency.  This  segment  then  uses  the  curve  selection 
routine,  the  parameter  values,  and  the  independent  variable  value 
to  determine  the  sum  of  square  deviations  of  the  data  being  con¬ 
sidered.  This  segment  then  returns  control  to  its  calling  routine. 

5.12  FUNCTION  FCT 

This  segment  is  used  to  select  the  curves  to  fit  the  data. 

The  parameter  values  and  independent  variable  value  are  used  to 
select  either  the  modulus  value  or  the  loss  factor  value  for  the 
parameter  values  and  independent  variable  value  being  considered 
This  segment  then  returns  a  value  to  its  calling  routine. 

5.13  FUNCTION  ALPHT 

This  segment  uses  the  temperature  of  the  data  point  along 
with  the  reference  temperature  for  the  material  to  calculate  the 
value  of  the  aT  curve  presented  in  Section  2.  This  aT  value  is 
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then  multiplied  by  the  frequency  of  the  data  point  to  determine 
the  independent  variable  value  for  the  data  point,  namely, 
reduced  frequency,  which  is  then  used  in  the  modulus  and  loss 
factor  curves  and  curve  fits.  This  segment  then  returns  the  aT 
value  to  its  calling  routine. 

5.14  FUNCTION  FI 

This  segment  of  the  program  is  used  to  evaluate  the 
modulus  curve  for  the  given  parameter  values  and  the  independent 
variable  value.  If  any  of  the  parameter  values  are  outside  the 
bounds  of  the  constraints  presented  in  Section  2,  the  value  of  the 
modulus  curve  is  set  to  a  large  number  in  hope  that  the  minimiza¬ 
tion  process  will  bring  the  parameter  values  within  the  constraint 
bounds  and  obtain  the  appropriate  curve  fit  to  the  modulus  data. 
This  segment  then  returns  the  modulus  curve  value  to  its  calling 
routine . 

5.15  FUNCTION  F2 

This  segment  of  the  program  is  used  to  evaluate  the  loss 
factor  curve  for  the  given  parameter  values  and  the  independent 
variable  value.  If  any  of  the  parameter  values  are  outside  the 
bounds  of  the  constraints  presented  in  Section  2,  the  value  of 
the  loss  factor  curve  is  set  to  a  large  number  in  the  hope  that  the 
minimization  process  will  bring  the  parameter  values  within  the 
constraint  bounds  and  obtain  the  appropriate  curve  fit  to  the 
loss  factor  data.  This  segment  then  returns  the  loss  factor  value 
to  its  calling  routine. 

5.16  FUNCTION  AL0G1 

This  segment  of  the  program  is  used  to  evaluate  the 
logarithm,  to  base  ten,  of  a  number.  The  HP-2100  mini-computer  has 
this  routine  available  as  a  utility,  but  has  been  shown  to  have 
several  bugs  in  the  routine.  Therefore,  to  ensure  proper 
operation,  this  short  segment  is  added  to  the  program.  This 
segment  then  returns  the  value  to  its  calling  routine. 
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5.17  SUBROUTINE  VERTL 

This  segment  of  the  program  is  used  to  output  a  string  of 
characters,  for  graph  labeling  purposes,  in  the  vertical  direc¬ 
tion.  The  graphics  package  needed  for  the  program  does  not  have 
a  routine  to  output  a  vertical  string  of  characters  on  the 
terminal.  This  segment  is  implemented  to  accomplish  this  by  using 
calls  to  routines  in  the  graphics  package.  This  segment  then 
returns  control  of  the  program  to  its  calling  routine. 

5.18  SUBROUTINE  SEARC 

This  segment  is  used  to  minimize  the  function  of  the 
parameters,  namely,  the  sum  of  square  deviations  or  least  squares. 
This  segment  is  a  package  routine  which  uses  the  direct  search 
technique  to  minimize  the  function.  This  segment  uses  the 
function  routine  to  evaluate  the  function  as  this  segment  searches 
directly  in  the  parameter  space  to  determine  the  parameter  values 
that  minimize  the  function  being  considered.  This  segment  then 
returns  control  to  its  calling  routine  after  this  segment  has 
determined  the  parameter  values  in  the  parameter  space  that 
minimizes  the  sum  of  square  deviations. 

5.19  SUBROUTINE  GAXES 

This  segment  of  the  program  generates  on  the  screen  of 
the  terminal  the  characteristic  graph  for  the  material  data. 
Through  extensive  calls  to  the  graphics  package  routines,  this 
segment  clears  the  screen,  displays  the  appropriate  grid  to 
accommodate  the  data,  graphs  the  modulus  data,  graphs  the  modulus 

curve,  graphs  the  loss  factor  data,  graphs  the  loss  factor  curve, 

2 

graphs  the  temperature  lines  ,  outputs  the  appropriate  scales  for 
each  of  the  axes,  and  labels  each  of  the  four  axes.  This  segment 
pauses  execution  of  the  program  so  the  user  can  examine  the  graph 
and  obtain  copies  of  the  graph  for  further  reference  and  returns 
control  of  the  program  to  its  calling  routine. 
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5.20 


SUBROUTINE  PLTPK 


w 


This  segment  of  the  program  is  used  to  output,  for 
display  by  the  user,  the  curves  to  which  the  data  was  fit  along 
with  the  value  of  the  reference  temperature  and  values  for  the 
parameters  of  the  curve  fits.  This  segment  also  outputs  the 
material  code  and  material  name  along  with  the  final  curve  and 
parameter  values  which  are  represented  by  the  characteristic 
graph.  This  segment  pauses  execution  of  the  program  for  the 
user  to  obtain  copies  of  the  screen,  and  then  returns  control  to 
its  calling  routine. 

5.21  SUMMARY 

Mentioned  in  the  above  subsections  are  only  brief 
descriptions  of  the  subprograms  that  form  the  complete  program. 
Any  detailed  description  would  require  a  line-for-line  analysis 
of  the  code  for  each  subprogram  to  thoroughly  examine  the 
entire  computer  program. 


i 


SECTION  VI 

USER  MANUAL  WITH  DESCRIPTIVE  SAMPLE  INPUT/OUTPUT 

This  section  presents  the  user  interaction  with  the  program 
to  generate  the  desired  representation  of  the  experimental  visco¬ 
elastic  material  data.  This  section  is  divided  into  two  subsections. 
The  first  subsection  gives  the  user  a  general  walk-through  of  the 
program  and  describes  the  commands  the  user  must  input  from  the 
terminal.  The  second  subsection  presents  sample  outputs  of  the 
program  along  with  a  descriptive  note  indicating  the  input  by  the 
user . 

The  user  is  to  note  that  all  underscored  items  which  follow 
indicate  user  input  to  the  program.  The  underscored  upper-case 
letters  must  be  entered  as  indicated  to  ensure  proper  format 
within  the  program.  The  lower-case  letters  indicate  free-format 
input  by  the  user  to  the  program.  In  this  case  the  items  may 
be  entered  as  desired  by  the  user.  If  several  underscored  items 
appear  in  the  user  manual  on  a  single  line,  the  user  only  needs 
to  separate  these  free-format  items  by  a  comma  or  a  space  to 
ensure  proper  entry  to  the  program.  After  each  input,  the  user 
enters  the  carriage  return-line  feed  sequence  (CR/LF)  to  enter 
the  input  to  the  program. 

6.1  USER  MANUAL 

The  user  must  first  initiate  execution  of  the  program.  To 
accomplish  this,  the  user  must  input: 

:RUN,  VIBR2  CR/LF  .  (1) 

The  program  begins  by  clearing  the  screen  and  outputs  the 
following : 

LIST  THE  STORED  DATA? (Y ,N)  (2) 

for  which  the  user  inputs  either: 

Yes  CR/LF  or  No  CR/LF  . 

If  No  is  entered  for  (2),  the  program  proceeds  to  Step  2. 

If  Yes  is  entered  for  (2) ,  the  program  outputs: 

ENTER  THE  MAXIMUM  NUMBER  OF  MATERIAL  CODES  (3) 
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for  which  the  user  inputs : 

max  CR/LF , 

the  maximum  number  of  material  codes  available  to  the  user 
(currently  200  3-14-78)  .  The  program  searches  the  data  files 

created  and  outputs  the  material  code,  material  name,  and  data 
sources  for  all  stored  material  data,  and  then  proceeds  to  Step  2. 


STEP  2 

The  program  responds: 

READ  IN  NEW  DATA? (Y ,N)  (4) 

for  which  the  user  inputs  either: 

Yes  CR/LF  or  No  CR/LF. 

If  No  is  entered  for  (4) ,  meaning  the  user  wants  to  retrieve 
data  already  stored  on  the  data  files,  the  program  outputs: 

ENTER  THE  CODE  OF  RETRIEVAL  MATERIAL  (5) 

for  which  the  user  inputs: 

m  CR/LF , 

where  m  is  the  material  code  for  the  material  data  to  be  retrieved. 
If  the  material  code  does  not  match  the  material  data  stored  on  the 
files,  the  program  outputs: 

MISMATCH  IN  RETRIEVED  MATERIAL  CODES  (6) 

and  proceeds  back  to  output  (5)  until  a  match  for  retrieved  data 
occurs  upon  which  the  program  proceeds  to  Step  3. 

If  Yes  is  entered  for  (4)  the  program  prepares  for  the  entry 
of  new  data  to  the  program. 

The  program  outputs : 

ENTER  HEADER  DATA 

ENTER  THE  CODE  OF  THE  MATERIAL  (7) 

for  which  the  user  inputs: 

i  CR/L^ , 

where  the  number  i,  between  1  and  200  inclusive,  corresponds  to  the 
position  the  following  data  is  to  be  stored  on  the  data  files. 
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The  program  then  outputs : 

ENTER  THE  MATERIAL  NAME  (MAX  60  CHAR.)  (8) 

for  which  the  user  inputs: 

name  CR/LF , 

the  name  of  desired  material. 

The  program  then  outputs: 

ENTER  THE  MANUFACTURER  DATA  SOURCE  (MAX  60  CHAR.) 

(9) 


for 

the 

The 

for 

the 

The 

for 

any 

The 


for 

the 

The 

for 


which  the  user  inputs: 

source  CR/LF, 

manufacturer  source  of  material  data, 
program  then  outputs: 

ENTER  THE  AFML  DATA  SOURCE  (MAX  60  CHAR.)  (10) 
which  the  user  inputs: 

source  CR/LF , 

AFML  source  of  material  data, 
program  then  outputs : 

ENTER  OTHER  DATA  SOURCES  (MAX  60  CHAR.)  (11) 

which  the  user  inputs : 

sources  CR/LF , 

other  sources  of  material  data, 
program  then  outputs : 

ENTER  THE  NUMBER  OF  DATA  POINTS  (MAX  64  PTS . ) 

(12) 


which  the  user  inputs : 

n  CR/LF , 

number  of  data  points  to  input, 
program  then  outputs : 

IS  THE  DATA  TO  BE  ENTERED  METRIC? (Y,N)  (13) 

which  the  user  inputs  either : 

Yes  CR/LF  or  No  CR/LF. 


If  Ye_s  is  entered  for  (13)  ,  a  flag  is  set  to  assume  the  data 
to  be  entered  is  metric. 
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If  No  is  entered  for  (13) ,  a  flag  is  set  to  assume  the 
data  to  be  entered  is  English  and  will  be  converted  to  metric 
immediately  after  entering. 

The  program  then  outputs: 

IS  THE  MODULUS  DATA  SHEAR? (Y,N)  (14) 

for  which  the  user  inputs  either: 

Yes  CR/LF  or  No  CR/LF. 

If  Yes  is  entered  for  (14) ,  a  flag  is  set  to  indicate  the 
data  is  for  shear  modulus  of  the  material. 

If  No  is  entered  for  (14)  ,  a  flag  is  set  to  indicate  the 
data  is  for  Young's  modulus  of  the  material. 

The  program  then  outputs  several  lines  for  the  user,  to 
indicate  the  program  is  now  ready  to  enter  the  viscoelastic 
material  data. 

The  user  inputs: 

i,  Mi,  n±»  tif  fL,  mi  CR/LF, 

th 

for  i  between  1  and  n  inclusive,  where  i  is  the  i  data  point, 
is  the  modulus,  ni  is  the  loss  factor,  t^  is  the  temperature, 
f^  is  the  frequency,  and  rtu  is  the  mode  number  for  that  data  point. 
The  program  converts  the  data  to  metric  if  entered  in  English  since 
all  computations  are  done  in  metric.  The  program  then  proceeds  to 
Step  3 . 

STEP  3 

The  program  clears  the  screen  and  proceeds  to  output  the 
new  or  retrieved  material  data  including  headings  and  labels,  where 
all  the  units  are  indicated  as  metric  since  all  computations  are 
done  in  metric.  The  program  then  pauses  so  that  the  user  can  obtain 
copies  of  the  screen  for  future  reference.  The  program  outputs: 
VIBR2 :  PAUSE  0000 

VIBR2  SUSP.  (15) 

@ 

to  continue  the  program  when  the  user  is  ready  the  user  inputs: 

: GO  CR/LF. 
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(16) 


r 


The  program  outputs : 

EDIT  THE  DATA? { Y  ,  N ) 

for  which  the  user  inputs  either: 

Yes  CR/LF  or  No  CR/LF . 

If  No  is  entered  for  (16),  the  program  proceeds  to  Step  4. 
If  Yes  is  entered  for  (16)  ,  the  program  proceeds  to  edit 
the  data  shown  on  the  screen.  The  editing  portion  of  the  program 
outputs : 

EDIT  THE  ABOVE  DATA:  ADD,  DELETE,  OR  END? (A , D , E) 

(17) 


for  which  the  user  inputs: 

Add  CR/LF  or  Delete  CR/LF  or  End  CR/LF. 

If  Add  is  entered  for  (17),  the  program  outputs: 

ENTER  THE  NUMBER  OF  DATA  POINTS  TO  ADD  (18) 

for  which  the  user  inputs: 

m  CR/LF , 

the  number  of  data  points  to  add  to  the  data  shown.  The  total  number 
of  data  points  cannot  exceed  64. 

The  program  then  outputs : 

IS  THE  DATA  TO  BE  ADDED  METRIC? (Y,N)  (19) 

for  which  the  user  inputs  either: 

Yes  CR/LF  or  No  CR/LF. 

If  Yes  is  entered  for  (19) ,  a  flag  is  set  to  indicate  the  data 
to  be  added  is  metric. 

If  No  is  entered  for  (19) ,  a  flag  is  set  to  indicate  the  data 
to  be  added  is  English  and  a  conversion  to  metric  will  be  performed 
after  the  data  is  entered. 


The  program  then  outputs: 

READY  TO  ADD  DATA  FOR  MATERIAL  IN  FORMAT  SHOWN  ABOVE 


for  which  the  user  inputs : 


i#  ni»  f^,  nu  CR/LF  , 

for  i,  between  1  and  m  inclusive,  being  the  ifc^  added  data  point 
for  the  material  data.  The  program  makes  the  conversion  to  metric 
if  necessary.  The  program  then  proceeds  to  output  (17)  again. 

If  Delete  is  entered  for  (17),  the  program  outputs: 

ENTER  THE  NUMBER  OF  DATA  POINTS 
TO  DELETE  AND  THE  INDICES  OF  THE  ROWS  TO  DELETE  (21) 

for  which  the  user  inputs : 

m,  i.  ,  i_,  ...  i  CR/LF  , 

_ l  i _ m 

and  the  specified  i_.  row  in  the  data  shown  on  the  screen  will  be 
deleted  from  the  list.  The  program  then  proceeds  to  output  (17) 
again. 

If  End  is  entered  for  (17) ,  the  program  terminates  the 
edit  sequence  and  proceeds  to  Step  4. 

STEP  4 

The  program  clears  the  screen,  produces  the  frequency  vs. 
temperature  plot  and  loss  factor  vs.  temperature  plot,  and  lists 
the  data  for  viewing  by  the  user.  A  computer  selected  Tq  is 
indicated  and  at  this  time  the  user  can  examine  the  data  closely. 

The  program  will  finally  output  (16)  again  and  the  user 
has  the  same  options  to  edit  the  data  as  in  Step  3.  If  the  user 
chooses  not  to  edit,  at  this  point  the  program  pauses  and  suspends 
by  outputting: 

VIBR2 :  PAUSE  0000 

VIBR2  SUSP  (22) 

0 
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at  this  time  the  user  can  obtain  copies  of  the  screen  for  future 
reference.  To  continue  the  program  the  user  inputs: 

:G0  CR/LF, 

and  proceeds  to  Step  5 . 

STEP  5 

To  begin  the  curve  fit  of  the  material  data  the  program 
outputs : 

CHANGE  TO  VALUE  OR  PARAMETER  VALUES? (Y,N)  (23) 

for  which  the  user  inputs: 

Yes  CR/LF  or  No  CR/LF. 

If  N<d  is  entered  for  (24),  the  program  proceeds  to  Step  6. 

If  Yes  is  entered  for  (24)  ,  the  program  outputs: 

ENTER  TO  VALUE  (24) 

for  which  the  user  inputs: 

t  CR/LF , 

the  t  value  for  the  material  to  the  nearest  10°C  by  either  user 
judgment  or  prior  experience  when  testing  the  material.  The  program 
then  proceeds  to  Step  6. 

STEP  6 

The  program  now  fits  the  material  modulus  data  to  a 
specified  curve.  The  program  clears  the  screen  and  outputs: 

FIT  MODULUS  DATA 

ENTER  THE  MINIMUM  VALUE  OF  FROM  PARAMETER  (25) 
for  which  the  user  inputs: 

x  CR/LF  , 

the  minimum  value  of  the  reduced  frequency  (FROM)  parameter  in  Hz. 
This  parameter  corresponds  to  the  value  of  F  T  at  the  inflection 
point  of  the  modulus  curve.  The  program  then  outputs: 

ENTER  THE  MINIMUM  VALUE  OF  MROM  PARAMETER  (26) 
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for  which  the  user  inputs : 

x  CR/LF, 

2 

the  minimum  value  of  the  modulus  (MROM)  parameter  in  N/m  .  This 
parameter  corresponds  to  the  value  of  the  modulus  curve  at  the 
inflection  point. 

The  program  then  outputs: 

ENTER  THE  MAXIMUM  VALUE  OF  MROM  PARAMETER  (27) 

for  which  the  user  inputs: 

x  CR/LF, 

o 

the  maximum  value  of  the  modulus  (MROM)  parameter  in  N/m  .  The 
program  then  outputs : 

ENTER  THE  MINIMUM  VALUE  OF  ML  PARAMETER  (28) 

for  which  the  user  inputs: 

x  CR/LF , 

2 

the  minimum  value  of  the  modulus  (ML)  parameter  in  N/m  .  This 
parameter  corresponds  to  the  lower  horizontal  asymptote  of  the 
modulus  curve. 

The  program  proceeds  to  fit  the  modulus  data  to  the  specified 
modulus  curve.  The  program  then  outputs  the  equation  of  the  modulus 
curve,  the  parameter  names,  and  the  values  of  the  parameters  that 
most  closely  describe  the  curve  fit  for  the  modulus  data.  The 
program  then  outputs: 

ACCEPT  THESE  VALUES? (Y,N)  (29) 

for  which  the  user  inputs  either: 

Yes  CR/LF  or  No  CR/LF. 

If  Yes  is  entered  for  (29) ,  the  program  proceeds  to  Step  8. 

If  No  is  entered  for  (29) ,  the  program  outputs: 

HOLD  ANY  PARAMETERS  CONSTANT  ?(Y,N)  (3Q) 

for  which  the  user  inputs: 

Yes  CR/LF  or  No  CR/LF. 
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If  No  is  entered  for  (30),  the  program  proceeds  to  Step  7. 

If  Yes  is  entered  for  (30),  the  program  outputs: 

ENTER  THE  NUMBER  AND  INDICES  OF  THE  CONSTANT 
PARAMETERS  (31) 

for  which  the  user  inputs: 

n,  •••  CR/LF t 

where  the  i^  are  the  indices  of  the  array  elements  that  correspond 
to  the  indicated  modulus  parameters.  Through  experience,  the  user 
may  want  to  hold  parameters  constant  for  a  given  material,  there¬ 
fore  this  option  is  available  to  the  user.  The  program  proceeds 
to  Step  7. 

STEP  7 

The  program  outputs: 

ENTER  FROM (Al)  PARAMETER  VALUE  (32) 

for  which  the  user  inputs: 

x  CR/LF  , 

the  user  defined  constant  or  variable  value  of  the  modulus  (FROM) 
parameter. 

The  program  then  outputs: 

ENTER  MROM (A2)  PARAMETER  VALUE  (33) 

for  which  the  user  inputs: 

x  CR/LF  , 

the  user  defined  constant  or  variable  value  of  the  modulus  (MROM) 
parameter.  The  program  then  outputs: 

ENTER  N(A3)  PARAMETER  VALUE  (34) 

for  which  the  user  inputs: 

x  CR/LF » 

the  user  defined  constant  or  variable  value  of  the  modulus  (N) 
parameter.  This  parameter  corresponds  to  the  slope  of  modulus  curve 
at  the  point  of  inflection. 
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The  program  then  outputs: 


ENTER  ML (A4 )  PARAMETER  VALUE  (35) 

for  which  the  user  inputs: 

X  CR/LF , 

the  user  defined  constant  or  variable  value  of  the  modulus  (ML) 
parameter.  The  program  proceeds  to  curve  fit  the  modulus  data  with 
these  user  defined  constant  and/or  variable  modulus  parameter  values. 
The  program  then  outputs  the  equation  of  the  modulus  curve,  the 
parameter  names,  and  the  values  of  the  parameters  that  now  most 
closely  describe  the  curve  fit  for  the  modulus  data.  The  program 
then  proceeds  to  output  (29)  again  for  the  user  option  input. 


STEP  8 

Next  the  program  fits  the  material  loss  factor  data  to 
a  specified  curve.  The  program  clears  the  screen  and  outputs: 

FIT  LOSS  FACTOR  DATA 

ENTER  THE  MINIMUM  VALUE  OF  ETAFROL  PARAMETER 

(36) 


for  which  the  user  inputs : 

x  CR/LF  , 

the  minimum  value  of  the  loss  factor  (ETAFROL)  parameter.  This 
parameter  corresponds  to  the  value  at  the  peak  of  the  loss  factor 
curve.  The  program  then  outputs: 

ENTER  THE  MINIMUM  VALUE  OF  SL  PARAMETER  (37) 

for  which  the  user  inputs: 

x  CR/LF  , 

the  minimum  value  of  the  slope  (SL)  parameter.  This  parameter 
corresponds  to  the  slope  of  the  asymptotic  line  for  the  low  values 
of  reduced  frequency  for  the  loss  factor  curve.  The  program  then 
outputs : 
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ENTER  THE  MINIMUM  VALUE  OF  FROL  PARAMETER 


(39) 


for  which  the  user  inputs: 

x  CR/LF , 

the  minimum  value  of  the  loss  factor  (FROL)  parameter  in  Hz. 

This  parameter  corresponds  to  the  reduced  frequency  value  at  the 
peak  of  the  loss  factor  curve. 

The  program  proceeds  to  fit  the  loss  factor  data  to  the 
specified  loss  factor  curve.  The  program  then  outputs  the  equation 
of  the  loss  factor  curve,  the  parameter  names,  and  the  values  of 
the  parameters  that  most  closely  describe  the  curve  fit  for  the 
loss  factor  data.  The  program  then  outputs: 

ACCEPT  THESE  VALUES? (Y,N)  (40) 

for  which  the  user  inputs  either: 

Yes  CR/LF  or  No  CR/LF  . 

If  Yes  is  entered  for  (40),  the  program  proceeds  to  Step  10. 
If  No  is  entered  for  (40),  the  program  outputs: 

HOLD  ANY  PARAMETERS  CONSTANT? (Y ,N)  (41) 

for  which  the  user  inputs  either: 

Yes  CR/LF  or  No  CR/LF. 

If  No  is  entered  for  (41) ,  the  program  proceeds  to  Step  9. 

If  Yes  is  entered  for  (41) ,  the  program  outputs: 

ENTER  THE  NUMBER  AND  INDICES  OF  THE  CONSTANT 
PARAMETERS  (42) 

for  which  the  user  inputs ; 

n,  i.^*  ij »  ...  CR/LF  , 

where  the  ij  are  the  indices  of  the  array  elements  that  correspond 
to  the  indicated  loss  factor  parameters.  Through  experience  the 
user  may  want  to  hold  parameters  constant  for  a  given  material, 
therefore  this  option  is  available  to  the  user.  The  program 
proceeds  to  Step  9. 
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STEP  9 


The  program  outputs : 

ENTER  ETAFROL (Bl)  PARAMETER  VALUE  (43) 

for  which  the  user  inputs: 

x  CR/LF  , 

the  user  defined  constant  or  variable  value  of  the  loss  factor 
(ETAFROL)  parameter.  The  program  then  outputs: 

ENTER  SL(B2)  PARAMETER  VALUE  (44) 

for  which  the  user  inputs: 

x  CR/LF  , 

the  user  defined  constant  or  variable  value  for  loss  factor  (SL) 
parameter.  The  program  then  outputs: 

ENTER  SH (B3 )  PARAMETER  VALUE  (45) 

for  which  the  user  inputs: 

x  CR/LF  , 

the  user  defined  constant  or  variable  value  for  the  loss  factor 
(SH)  parameter.  The  program  then  outputs: 

ENTER  FROL (B4)  PARAMETER  VALUE  (46) 

for  which  the  user  inputs : 

x  CR/LF  , 

the  user  defined  constant  or  variable  value  of  the  loss  factor 
(FROL)  parameter.  The  program  then  outputs: 

ENTER  C(B5)  PARAMETER  VALUE  (47) 

for  which  the  user  inputs : 

x  CR/LF  , 


the  user  defined  constant  or  variable  value  of  the  loss  factor  (C) 


parameter.  This  parameter  corresponds  to  the  shape  of  the  peak  of 
the  loss  factor  curve  (larger  C  implies  flatter  peak,  smaller  C 
implies  sharper  peak) . 
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The  program  proceeds  to  curve  fit  the  loss  factor  data  with  these 
user  defined  constant  and/or  variable  loss  factor  parameter  values. 
The  program  then  outputs  the  equation  of  the  loss  factor  curve, 
the  parameter  names,  and  the  value  of  the  parameters  that  now  most 
closely  describe  the  curve  fit  for  the  loss  factor  data.  The 
program  then  proceeds  to  output  (40)  again  for  the  user  option 
input. 

STEP  10 


The  program  outputs: 

ENTER  THE  INTERVAL  OF  TEMPERATURE 

LINES  (48) 

for  which  the  user  inputs: 

x  CR/LF 

the  interval  of  the  displayed  temperature  lines.  If  the  interval 
value  entered  is  less  than  20°C,  the  half-band  width  is  assigned 
to  100°C.  If  the  interval  value  entered  is  greater  than  20°C, 
the  half-band  width  is  assigned  to  200°C. 

The  desired  graphs  are  then  generated.  Since  the  program 
has  determined  the  values  of  the  parameters  for  each  curve  and  the 
user  has  accepted  these  values,  the  program  clears  the  screen  and 
generates  the  desired  grid.  The  grid  is  generated  in  such  a  way 
that  all  the  data  points  will  fall  within  its  bounds.  The  data 
points  for  both  the  modulus  and  loss  factor  are  plotted  as  a 
function  of  reduced  frequency.  The  curves  for  both  are  plotted 
using  the  parameters  determined  and  being  a  function  of  reduced 
frequency.  The  axes  are  labeled  with  the  correct  label  in  metric 
units  along  with  the  scale  for  each  axis. 

The  final  plot  is  now  complete  and  the  program  pauses  and 
suspends  by  outputting: 


VIBR2:  PAUSE  0000 
VIBR2  SUSP 
@ 
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(49) 


the  user  can  obtain  copies  of  the  screen  for  future  reference .  To 
continue  the  program  the  user  then  inputs : 


;G0  CR/LF  . 

The  program  clears  the  screen  and  outputs  the  material  name, 
equations,  parameter  names,  parameter  values,  and  proceeds  to  pause 
and  suspend  by  outputting: 

VIBR2:  PAUSE  0000 

VIBR2  SUSP  (50) 

e 

the  user  can  obtain  copies  of  the  screen  which  identify  the  previous 
plot.  To  continue  the  program  the  user  then  inputs: 

:GO  CR/LF  . 


STEP  11 

The  program  outputs  (23)  again  and  the  user  has  the  options 
of  changing  t  and/or  the  parameters  given  to  define  the  relationship 
of  the  curves  generated  for  the  experimental  data. 

If  Yes  is  entered  for  (23),  the  program  goes  to  output  (24) 
and  proceeds  from  that  point. 

If  No  is  entered  for  (23),  the  program  outputs: 

STORE  THE  DATA? (Y ,N)  (51) 

for  which  the  user  inputs  either: 

Yes  CR/LF  or  No  CR/LF  . 

If  No  is  entered  for  (51),  the  program  proceeds  to  Step  12. 

If  Yes  is  entered  for  (51),  the  program  outputs: 

ABOVE  PLOTTED  DATA  AND 

PARAMETERS  ARE  STORED  ON  (52) 

FILES:  VIBR3,  VIBR4 ,  AND  VIBR5 

to  indicate  to  the  user  the  data  has  been  stored  on  the  designated 
files  for  further  reference  at  another  time.  The  program  then 
pauses  and  suspends  by  outputting: 
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VIBR2:  PAUSE  0000 

VIBR2  SUSP 


@  (53) 

for  which  the  user  inputs: 

:GO  CR/LF, 

to  continue  the  program  to  Step  12. 

STEP  12 

The  program  outputs: 

BEGIN  ANOTHER  PROBLEM? (Y,N)  (54) 

for  which  the  user  inputs  either: 

Yes  CR/LF  or  No  CR/LF. 

If  Yes  is  entered  for  (54) ,  the  program  proceeds  to  output 
(2)  and  the  user  begins  another  problem  (i.e.,  enter  or  retrieve 
another  material  data  set,  etc.). 

If  No  is  entered  for  (54),  the  program  terminates. 

6.2  Sample  Input/Output 

The  following  pages  contain  example  input/output  for  the 
program  user  manual  described  above.  All  the  options  cannot  bo 
presented  in  that  the  number  of  pages  would  be  tremendous,  so 
just  brief  examples  of  the  input/output  are  presented  to  illustrate 
the  features  present  in  the  program  and  described  in  the  above  user 
manual . 
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SECTION  VII 


SUMMARY  AND  RECOMMENDATIONS 

The  computer  program  described  in  this  report  is  operational 
on  the  computer  equipment  at  AFML.  The  program  is  being  used  to 
reduce,  edit,  and  graphically  represent  viscoelastic  material 
properties  data  for  further  development  of  a  standardized  material 
properties  handbook. 

The  program  currently  requires  practically  the  entire 
available  core  of  the  computer  equipment  when  executing  for  the 
user.  Therefore,  any  modifications  to  the  computer  program  may 
be  costly  if  they  require  the  use  of  more  core  than  is  currently 
available  for  the  equipment.  If  modifications  are  required 
which  would  increase  the  program  size  beyond  the  available  core, 
the  program  would  have  to  be  split  into  suitable  segments.  By 
splitting  the  program  into  suitable  segments,  it  could  be  executed 
in  the  available  core  of  the  equipment.  Segmentation  allows  for 
the  execution  of  segments  as  needed  by  the  program,  thus  reducing 
the  overall  size  of  the  executing  program.  Segmentation  is  avail¬ 
able  on  the  AFML  equipment. 

Future  recommendations  for  this  computer  program  would  be 
to  give  the  user  the  option  of  defining  or  using  one  of  several 
curve-fit  equations  for  the  modulus  and/or  loss  factor  data. 

On  this  case,  the  user  could  select  the  equations  that  he  feels 
best  describe  the  material  data.  Also,  the  user  should  have  the 
option  of  choosing  a  scheme  to  find  the  reference  temperature  (tQ) 
for  the  material,  and  in  gome  way  to  find  the  shift  parameter  (at) 
curve  that  best  describes  the  shift  needed  for  each  material. 

These  capabilities  if  added  to  the  program  could  give  the  user 
a  better  representation  of  the  viscoelastic  material  properties 
data. 

Several  further  considerations  are  the  implementation  of 
features  into  this  program,  or  possible  other  programs,  to  use  the 
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computer  to  control  the  experimental  test  facility.  The  use  of  the 
computer  for  all  data  recording,  data  processing,  data  editing,  and 
data  representation  would  reduce  the  time  needed  to  collect  the 
data,  and  also  take  over  many  more  manual  operations  currently  being 
performed.  Through  this  major  task,  reported  herein,  the  use  of 
the  computer  for  all  the  data  processing  operations  has  become  more 
feasible.  Future  developments  include  the  integration  of  the  com¬ 
puter  in  to  the  test  facility,  and  programs  implemented  on  the  com¬ 
puter  to  assist  the  user  to  effectively  record  data  from  the  test 
facility,  and  proceed  through  the  representation  of  the  data  by  the 
computer.  The  amount  of  time  spent  by  the  user,  therefore,  would 
be  reduced,  because  the  user  now  will  have  a  fast  digital  computer 
to  accomplish  in  minutes  many  of  the  operations  the  user  would 
currently  do  in  several  hours. 
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APPENDIX  A 

PARAMETER  NAME  DESCRIPTIONS 


Following  is  a  list  of  the  parameter  names  used  for  the 
modulus  and  loss  factor  curves,  indicating  the  units,  and 
describing  the  parameter  for  the  appropriate  curve. 


Name 

Units 

TO 

°C 

Reference  temperature  to  shift  and  minimize 

scatter  for  modulus  and  loss  factor  curves 

<V* 

F 

HZ 

Frequency  of  observation  (f). 

FR 

HZ 

Frequency  of  observation  reduced  by  the  shift 
factor  induced  by  the  reference  temperature 
(fr  =  aTf) . 

T 

°C 

Temperature  of  observation  (t) . 

MODULUS 

CURVE : 

M 

N/m2 

Observed  modulus  for  material  (M) . 

FROM 

HZ 

Reduced  frequency  at  the  inflection  point  of 

the  modulus  curve  (f  _ ) . 

rom 

MROM 

N/m2 

Value  of  the  modulus  curve  at  the  inflection 

P°int  (Mrom) * 

N 

- 

Shape  parameter  for  the  modulus  curve  (N) . 

ML 

N/m2 

Value  of  the  lower  asymptote  for  the  modulus 
curve  (M^) . 

LOSS  FACTOR  CURVE: 

ETA  -  Observed  loss  factor  for  material  { n ) . 

ETAFROL  -  Value  of  the  loss  factor  curve  at  the  peak 

(nfroJl)  * 

Slope  of  asymptotic  line  for  low  values  of 
reduced  frequency  (S^) . 


SL 


Name 


Units 


SH 

FROL  HZ 

C 


Slope  of  asymptotic  line  for  high  values  of 
reduced  frequency  (S^) . 

Reduced  frequency  at  the  peak  of  the  loss 


factor  curve  (f  „) 
roZ 


Shape  parameter  for  the  loss  factor  curve 


(C)  . 


APPENDIX  B 
PROGRAM  LISTING 


This  section  contains  a  compiled  listing  of  the  program 
described.  This  program  has  been  compiled  on  the  CDC  computer 
system  at  Wright-Patterson  Air  Force  Base,  because  of  the  inability 
of  the  HP-2100  and  equipment  available  to  present  a  neat  and  clean 
compiled  listing  of  the  program.  There  are  no  differences  in  the 
FORTRAN  source  presented  herein  and  that  stored  on  the  HP-2100 
computer  system  with  permanent  file  name  VIBRl. 
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FTN4 


PROGRAM  VIBR2 

COMMON  DUMMY (60) , IR, IW, IBUF (128 ) ,NPTS , XIN ( 64 , 5) ,YD(64) , 

1  COFTO ( 3 ) , TEMPO ,CHISQ , NPAR , DELTA ( 5 ) , DELMI(  5) , 

2  MASK{ 5) ,GPAR( 4) ,ETPAR( 5) ,CNSTR( 4) 

C 

C  VIBR3=#MATERIALS  (BLOCKS) 

C  VIBR4=5*VIBR3 

C  VIBR5=VIBR3 

C 

I Y=1HY 

IR=1 

IW=1 

CALL  INITT( 1200 ) 

CALL  TERM (2 ,1024) 

50  CALL  ERAS 

CALL  HOME 
CALL  CHSIZ ( 3 ) 

CALL  ANMOD 
CALL  REDIN 
WRITE (I W, 1000) 

READ(I  R, 1005)  I YN 
IF(IYN.EQ.IY)  CALL  EDITO 
100  CALL  FITT0 

CALL  PLTT0 
WRITE(I W, 1000 ) 

READ(I  R, 1005 )  I YN 
IF(IYN.NE.IY)  GO  TO  300 
CALL  EDITO 
GO  TO  100 
300  PAUSE 

WRITE (I W, 1500 ) 

READ(I R, 1005 )  I  YN 
IF(IYN.EQ.IY)  GO  TO  450 
400  CALL  MODUS 

CALL  ETA 
WRITE(I W, 4000) 

READ (I R, * )  DELT0 
CALL  GAXES(DELTO) 

CALL  PLTPK ( 3 ) 

WRITE (I W, 1500 ) 

READ(I  R, 1005 )  I YN 
IF(IYN.NE.IY)  GO  TO  500 
450  WRITE (I W, 1600 ) 

READ(I  R,  *  )  TEMPO 
GO  TO  400 

500  WRITE  (I  W,  3000  ) 

READ(IR,  1005)  I  YN 
IF(IYN.EQ.IY)  CALL  STORE 
WRITE(I  W, 2000 ) 

READ(I R, 1005 )  I YN 
IF(IYN.EQ.IY)  GO  TO  50 
CALL  FINIT( 0,100) 

1000  FORMAT( 2X , 19HEDIT  THE  DATA?(Y ,N) ) 

1005  FORMAT (1A1) 


78 


1500 

1600 

2000 

3000 

4000 


80 

90 


FORMAT (2X, 4 1HCHANGE  TO  VALUE  OR  PARAMETER  VALUES? (Y ,N ) ) 
FORMAT( 2X , 14HENTER  TO  VALUE) 

FORMAT ( 2X, 27HBEGIN  ANOTHER  PROBLEM? (Y ,N) ) 

FORMAT ( 2X , 20HSTORE  THE  DATA? (Y  ,N )  ) 

FORMAT (2X,21HENTER  THE  INTERVAL  OF 
1  /2X , 17HTEMPERATURE  LINES) 

END 


SUBROUTINE  REDIN 

COMMON  DUMMY (60) , IR, IW, IBUF ( 128 ) , NPTS , XIN( 64 , 5 ) ,YD(64) , 

1  COFT0( 3) , TEMPO ,CHISQ,NPAR, DELTA (5) ,DELMI(5) , 

2  MASK( 5) , GPAR( 4 ) ,ETPAR(5) ,CNSTR(4) 

DIMENSION  BUF(320) ,NAME( 3) 

IY=1HY 

WRITE(I W, 5000 ) 

READ(I R, 1005 )  I YN 
IF(IYN.NE.IY)  GO  TO  90 
WRITE (I W, 6000 ) 

READ  (I  R,  * )  NMAX 
ICR=14 
LIBUF=128 
ICON=3B 
NAME( 1 ) =2HVI 
NAME ( 2 ) =2HBR 
NAME( 3)=2H3 
DO  80  I®1 /NMAX 
INS=I -1 

CALL  EXEC (ICR, ICON , IBUF , LIBUF , NAME , INS ) 

IF(IBUF(121)  . NE .  I )  GO  TO  80 

WRITE  (I  W,  3200  )  IBUF  (121)  ,  (I  BUF  ( J  )  ,J=1,120) 

CONTINUE 
WRITE(I W, 1000 ) 

READ(I  R, 1005 )  I YN 

IF (I YN . NE . IY )  GO  TO  200 

WRITE (I  W, 2000 ) 

WRITE (I W, 2100 ) 

READ (I R, *)  MATL 
IBUF ( 121 )=MATL 
WRITE (I W, 2200 ) 

READ(I R, 1010 )  (I BUF (I )  ,1=1,30) 

WRITE(I W, 2300 ) 

READ(I R, 1010 )  (I BUF(I ) ,1=31,60) 

WRITE (I W, 2400 ) 

READ(I  R,  1010  )  (IBUF(I  )  ,1=61,90) 

WRITE (I W, 2500 ) 

READ(I  R,  1010  )  (IBUF(I  )  ,1=91,120) 

WRITE (I W, 2600 ) 

READ  (I  R,  *  )  NPTS 
IBUF(122)=NPTS 
METRC=0 
WRITE(I W, 2650) 
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READ(I R/ 1005 )  IYN 
IF(IYN.EQ.IY)  METRC® 1 
IBUF ( 124 ) =1 
WRITE (I W, 2675 ) 

READ(I R/ 1005 )  IYN 
IF  (IYN.EQ.IY)  I  BUF (124 )  =0 
WRITE (I W, 2700 ) 

WRITE (I W, 2800 ) 

IF (METRC. EQ. 1)  WRITE(I W, 2850 ) 

IF (METRO. EQ. 0)  WRITE(I W, 2900) 

DO  100  1=1 , NPTS 

100  READ(IR,*)  N , (XIN(I , J) , J=1 , 5) 

I F ( METRC . EQ . 1 )  GO  TO  500 
DO  150  1=1, NPTS 
XIN(I , 1) =XIN(I , 1) *  6894 . 757 
150  XIN (I ,3) =5. *(XIN(I , 3 ) -32. )/9. 

GO  TO  500 

200  WRITE  (I  W, 3000) 

READ(I  R, * )  MATL 
ICR=14 
LBUF=640 
LIBUF=128 
ICON=3B 
INS=MATL— 1 
IRS=5*I  NS 
NAME(1)=2HVI 
NAME( 2 ) =2HBR 
NAME( 3) =2H3 

CALL  EXEC (I CR , ICON ,IBUF,LIBUF, NAME , INS ) 

IF(IBUF(121) . EQ.MATL)  GO  TO  250 
WRITE (I W, 3100 ) 

GO  TO  200 
250  NAME(3)=2H4 

NPTS=IBUF(122) 

CALL  EXEC  (I  CR, ICON , BUF , LBUF , NAME , IRS ) 

DO  300  1=1, NPTS 
DO  300  J=1 , 5 

300  XIN  (I  ,  J )  =BUF  ( 5*  (I -1 ) +J ) 

500  CALL  ERAS 
CALL  HOME 
CALL  ANMOD 

WRITE ( I W, 3200 )  I BUF (121) , (IBUF(J) ,J=1,120) 

WRITE (I W, 2800) 

WRITE ( I W, 2850 ) 

DO  600  1=1, NPTS 

600  WRITE(I W, 3400 )  I , (XIN (I ,J) , J=1 , 5) 

PAUSE 

RETURN 

1000  FORMAT ( 2X , 22HREAD  IN  NEW  DATA? (Y ,N) ) 

1005  FORMAT (1A1) 

1010  FORMAT( 30A2) 

2000  FORMAT ( 2X , 17HENTER  HEADER  DATA) 

2100  FORMAT ( 2X , 30HENTER  THE  CODE  OF  THE  MATERIAL) 

2200  FORMAT( 2X , 37HENTER  THE  MATERIAL  NAME ( MAX  60  CHAR.)) 

2300  FORMAT (2X,48HENTER  THE  MANUFACTURER  DATA  SOURCE(MAX  60  CHAR.)) 
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2400  FORMAT( 2X , 40HENTER  THE  AFML  DATA  SOURCE(MAX  60  CHAR.)) 

2500  FORMAT ( 2X, 38HENTER  OTHER  DATA  SOURCES(MAX  60  CHAR.)) 

2600  FORMAT ( 2X, 44HENTER  THE  NUMBER  OF  DATA  POINTS(MAX  64  PTS.)) 
2650  FORMAT( 2X , 38HIS  THE  DATA  TO  BE  ENTERED  METRIC? (Y ,N ) ) 

2675  FORMAT ( 2X , 31HIS  THE  MODULUS  DATA  SHEAR? (Y ,N)  ) 

2700  FORMAT { 2X , 32HREADY  TO  ENTER  DATA  FOR  MATERIAL) 

2800  FORMAT ( 2X , 3HNO . , 3X , 7HMODULUS , 5X , 4HLOSS , 5X , 5HTEMP . ,4X,5HFREQ. , 
1  3X , 4HMODE ) 

2850  FORMAT ( 9X , 6HN/M**2 , 4X , 6HF ACTOR / 3X , 6HDEG .  C,5X, 

1  2HHZ , 5X , 3HNO . ) 

2900  FORMAT ( 7X , 8HLB A N**  2 , 4X , 6HFACTOR , 3X , 6HDEG .  F,5X, 

1  2HHZ ,  5X , 3HNO. ) 

3000  FORMAT( 2X , 36HENTER  THE  CODE  OF  RETRIEVAL  MATERIAL) 

3100  FORMAT { 2X , 36HMISMATCH  IN  RETRIEVED  MATERIAL  CODES) 

3200  FORMAT( 2X , 15HMATERIAL  CODE  :,I3, 

1  /2X , 10HMATERIAL  ., 30A2 ,/12X , 12HDATA  SOURCES , 

2  /2X , 14HMANUFACTURER  : , 30A2 , /2X , 6HAFML  :,30A2, 

3  /2X , 7HOTHER  :,30A2/) 

3400  FORMAT (2X,I3,1PE13. 5,0PF9. 4,F9. 1 ,  F9 . 1, F5. 0) 

5000  FORMAT( 2X, 26HLIST  THE  STORED  DATA? (Y ,N) ) 

6000  FORMAT( 2X , 42HENTER  THE  MAXIMUM  NUMBER  OF  MATERIAL  CODES) 

END 


SUBROUTINE  FITT0 

COMMON  DUMMY ( 60 ) , IR, IW, IBUF ( 128 ) ,NPTS , XIN ( 64 , 5 ) ,YD(64) , 

1  COFT0(3 ) , TEMPO, CHISQ,NPAR, DELTA (5) ,DELMI(5) , 

2  MASK ( 5 ) , GPAR( 4 ) , ETPAR( 5 ) ,CNSTR( 4 ) 

DIMENSION  X( 64 ) ,COF ( 3 ) 

DO  100  1=1 , NPTS 
X(I  )  =XIN (I  ,  3) 

ETA=XIN(I ,2) 

100  YD(I )=ALOGl(ETA) 

CALL  POLYF (X , COF ) 

DO  200  1=1,3 
200  COFTO  (I  )  =COF  (I  ) 

T0=-COF (2 )/ ( 2 . *  COF ( 3 ) ) 

IT0=T0/10.+. 5* ABS(T0 )/ TO 
TEMP0=IT0*10. 

RETURN 

END 


SUBROUTINE  POLYF (X, COF) 

COMMON  DUMMY ( 60 ) , IR, IW, IBUF {128 ) , NPTS, XIN( 64,5) ,YD(64) , 

1  COFTO ( 3 ) , TEMPO, CHISQ,NPAR,DELTA(5) ,DELMI(5) , 

2  MASK( 5 ) ,GPAR( 4 ) ,ETPAR(5) ,CNSTR(4) 

DIMENSION  X{ 64 ) ,COF ( 3 ) , A{ 3 , 3 ) ,B(3,3),C(3) 

DO  100  1=1,3 
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100  A(I,J)=0. 


A(1,1)=NPTS 
DO  200  1=1 ,NPTS 
A( 1 , 2 ) =A ( 1 , 2) +X(I  ) 

A  ( 1 , 3 )  =A  { 1 , 3 )  +X  (I  )**2 
A(2,3)=A(2,3)+X(I )**3 
A  ( 3 , 3 )  =A  (  3 , 3 )  +X  (I  )  *  *4 
C(1)=C  (1)+YD(I  ) 

C  (2 )  =C  ( 2 )  +X(I  )*YD(I  ) 

200  C  (  3 )  =C  (  3 ) +YD(I  )*X(I  )**2 

A(2,1)=A(1,2) 

A(2,2)=A(1,3) 

A(3,1)=A{1,3) 

A(3,2)=A(2,3) 

B(1,1)=A(2,2)*A(3,3)-A(3,2)*A(2,3) 
B(2,1)=A(3,1)*A<2,3)-A(2,1)*A<3,3> 
B(3,1)=A(2,1)*A(3,2) -A ( 3 , 1 )  *  A(  2 , 2 ) 
B(1,2)»A(3,2)*A(1,3)-A<1,2)*A(3#3) 
B(2/2)=A(l,l)*A(3,3)-A(3,l)*A(l/3) 
B(3,2)=A(l,2)*A(3,l)-A(l,l)*A(3/2) 
B(l,3)=A(l,2)*A{2/3)  -A  (2,2)*A(1,3) 
B(2,3)=A(2,l)*A(lf3)-A(l,l)*A(2,3) 
B(3,3)=A(1,1)*A(2,2)-A(1,2)*A(2,1) 

DET=A (1,1)*B(1,1)+A(1,2)*B(2,1) +A (1,3)*B(3,1) 
DO  300  1=1,3 
DO  300  J=1 , 3 

300  A  (I  ,  J )  =B  (I  ,  J ) /D  ET 
DO  400  1=1,3 
COF(I  )  =0 . 

DO  400  J=1 , 3 

400  COF(I  )=A(I  ,  J )  *  C(  J  ) +COF  (I  ) 

RETURN 

END 


SUBROUTINE  PLTT0 

COMMON  DUMMY ( 60 ) ,IR, IW, IBUF(128 ) ,NPTS , XIN( 64 , 5) ,YD(64) , 

1  COFTO ( 3 ) , TEMPO , CHI SQ , NPAR , DELTA ( 5 ) , DELMI  (  5 ) , 

2  MASK( 5) ,GPAR( 4 ) ,ETPAR(5) ,CNSTR(4) 

CALL  ERAS 

ETMAX=XIN( 1,2) 

ETMIN=XIN (1,2) 

TMAX=XIN (1,3) 

TMIN=XIN (1,3) 

FMAX=XIN{ 1,4) 

FMIN=XIN( 1,4) 

DO  10  1=1 ,NPTS 

IF ( XIN(I  ,2) .GT.ETMAX)  ETMAX=XIN(I  ,2) 

IF (XIN (I ,2) .LT.ETMIN)  ETMIN=XIN(I ,2) 

IF(XIN(I  ,  3  )  .  GT . TMAX )  TMAX=XIN(I  ,3) 
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IP (XIN(I  ,  3)  .  LT.TMIN)  TMIN«XIN(I  ,3) 

IF  (XIN(I  ,  4 )  . GT .  FMAX )  FMAX«XIN(I  ,4) 

IF (XIN(I ,4) . LT.FMIN)  FMIN«XIN(I ,4) 

10  CONTINUE 

I EXPN-ALOG1 (ETMIN ) 

ETMIN-IEXPN-1 
I EXPN-ALOG1 ( ETMAX ) 

ETMAX® I EXPN + 1 
IMULT-TMIN/10. 

TMIN*10. *(I MULT-1) 

IMULT*T MAX/10. 

TMAX*10 . * (I MULT+1 ) 

IEXPN=ALOGl (FMIN) 

FMIN-IEXPN-1 

IEXPN=ALOGl(FMAX) 

FMAX*IEXPN+1 

CALL  PLTT1 ( 4 , TMIN , TMAX ,  FMIN , FMAX ) 

CALL  PLTTl ( 2 , TMIN , TMAX , ETMIN , ETMAX ) 

CALL  RESET 
CALL  HOME 
CALL  CHSIZ ( 3 ) 

CALL  ANMOD 

WRITE(I  W,  3000)  I  BUF  ( 121 )  ,  (I  BUF(J)  , J=l, 30) 

WRITE(I W, 2800) 

WRITEd  W,  2850 ) 

DO  300  1=1 ,NPTS 

300  WRITE  (I  Wf  4000  )  I  ,  ( XIN  (I  ,  J )  ,  J=1 , 5) 

RETURN 

3000  FORMAT ( 2X , 15HMATERI AL  CODE  : , 13 ,/2X , 10HMATERIAL  :,30A2) 

2800  FORMAT ( 2X, 3HNO. , 3X , 7HMODULUS , 5X , 4HLOSS , 5X , 5HTEMP . ,4X,5HFREQ. , 
1  3X, 4HMODE) 

2850  FORMAT ( 9X , 6HN/M**2 , 4X , 6HFACTOR, 3X , 6HDEG .  C,5X, 

1  2HH  Z , 5X , 3HNO .  ) 

4000  FORMAT ( 2X,I3,1PE13. 5f0PF9. 4,F9. 1,F9. 1,F5.0) 

END 


SUBROUTINE  PLTTl  (N  UM , TMIN , TMAX , ETMIN , ETMAX ) 

COMMON  DUMMY ( 60 ) , IR, IW, IBUF ( 128 ) ,NPTS , XIN( 64 , 5) ,YD(64) , 

1  COFTO ( 3 ) , TEMPO, CHISQ,NPAR, DELTA (5) ,DELMI(5) , 

2  MASK ( 5 ) ,GPAR( 4 ) ,ETPAR(5) ,CNSTR(4) 

DIMENSION  11(2) , LABEL ( 9) 

ITEN=2H10 

CALL  DWIND(T MIN , TMAX , ETMIN , ETMAX ) 

IF (NUM . EQ. 4 )  CALL  TWIND( 560 , 1010 , 420 , 770 ) 

IF(NUM. EQ. 2)  CALL  TWIND( 560 , 1010 , 30 , 380 ) 

CALL  MOVEA  (TMIN,  ETMIN) 

CALL  DRAWA (TMAX, ETMIN) 

CALL  MOVEA (TMIN, ETMIN) 

CALL  DRAWA (TMIN, ETMAX) 

TM=TMIN 

80  CALL  MOVEA (TM, ETMIN) 
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CALL  DREL( 0 , 5) 

TM=TM+10. 

IF(TM.LE.TMAX)  GO  TO  80 
ETM=ETMIN 

90  CALL  MOVEA  (TMIN  ,  ETM ) 

CALL  DREL( 5,0) 

ETM=ETM+1 . 

IF (ETM . LE . ETMAX )  GO  TO  90 
CALL  CHSIZ ( 4 ) 

DO  100  1=1 ,NPTS 
TM=XIN (I ,3) 

ETM=XIN(I  ,  NUM) 

ETE=ALOGl(ETM) 

CALL  PNTA  (TMf  ETE) 

CALL  CODE 
WRITE( J , 1000 )  I 
CALL  MREL(2,-5) 

CALL  AOPST (  2 ,  J ) 

100  CONTINUE 

I F  (N  UM .  EQ .  4 )  GO  TO  115 
TM=TMIN 

ETM=( CDFT0 ( 1 ) +TM* ( COFTO ( 2 ) +TM*  COFTO ( 3 ) ) ) 
CALL  MOVEA  (TM, ETM) 

110  TM=TM+. 5 

ETM=(  COFTO  ( 1 )  +TM  *  (  COFTO  (  2  )  +TM  *  COFTO  (  3 )  )  ) 
CALL  DRAW  A  (TM,  ETM) 

IF(TM.LT.TMAX)  GO  TO  110 
115  CALL  CHSIZ (3) 

TM=TMIN 

MULT=(TMAX-TMIN)/100.+1. 

120  IT=TM 

CALL  CODE 
WRITE(I  1,2000)  IT 
CALL  MOVEA (TM,ETMIN) 

CALL  MREL( -24 ,-14 ) 

CALL  AOPST (4, II) 

TM=TM+MULT*10 . 

IF(TM.LE.TMAX)  GO  TO  120 
ETM=ETMIN 
130  IETME=ETM 
CALL  CODE 

WRITE  (J  ,1000)  IETME 
CALL  MOVEA  (TMIN,  ETM) 

CALL  MREL( -36 , -7 ) 

CALL  CHSIZ (3) 

CALL  AOPST ( 2 , ITEN ) 

CALL  MREL( 1,7) 

CALL  CHSIZ (4) 

CALL  AOPST ( 2, J) 

ETM=ETM+1 . 

IF(ETM.LE. ETMAX)  GO  TO  130 
CALL  CHSIZ (3) 

IF(NUM. EQ . 4 )  CALL  MABS( 754 , 390 ) 

IF (N UM. EQ . 2 )  CALL  MABS(754,0) 

LABEL ( 1 ) =2HTE 
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LABEL ( 2 ) =2HMP 
LABEL( 3 ) =2HER 
LABEL { 4 } =2HAT 
LABEL( 5 ) =  2HUR 
LABEL(6)=2HE 
LABEL( 7 ) =2HDE 
LABEL ( 8 ) =2HG. 

LABEL ( 9 ) =2H  C 

CALL  AOPST( 18, LABEL) 

IF(NUM.EQ.  4)  GO  TO  140 

CALL  MABS( 570 , 360 ) 

LABEL( 1)=2HET 

LABEL( 2 ) =1HA 

CALL  AOPST ( 3 , LABEL ) 

CALL  MABS( 900 , 360) 

LABEL( 1 ) =2HT0 

LABEL ( 2 ) =2H  = 

CALL  AOPST ( 4 , LABEL ) 
IT=T EMPO 
CALL  CODE 
WRITE  (1 1,2000)  IT 
CALL  AOPST (4, II) 

RETURN 

140  CALL  MABS ( 570 , 750 ) 
LABEL ( 1 ) =2HFR 
LABEL( 2)=2HEQ 
LABEL ( 3)=1H. 

CALL  AOPST (5, LABEL) 
RETURN 

1000  FORMAT (I  2 ) 

2000  FORMAT (14) 

END 


SUBROUTINE  EDIT0 

COMMON  DUMMY ( 60 ) , IR, IW, IBUF ( 128 ) ,NPTS , XIN(  64 , 5)  fYD(64) , 

1  COFT0( 3) , TEMPO, CHISQ,NPAR, DELTA (5) ,DELMI(5) , 

2  MASK( 5) , GPAR( 4 ) ,ETPAR(5) ,CNSTR(4) 

DIMENSION  IDE( 64 ) 

ID=1HD 

IE=1HE 

IY*1HY 

10  WRITE(IW,  1000) 

READ (I R, 1005 )  I AD 
IF (I AD. EQ« ID)  GO  TO  100 
IF  (I  AD.  EQ.  IE)  RETURN 
WRITE(IW,  1100) 

READ(I  R,  * )  NUM 
WRITE (IW,1150) 

READ(I R, 1005 )  IYN 
METRC-0 

IF  (IYN.  EQ.  IY)  METRC*1 
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NL=NPTS+1 
NH=NPTS+NUM 
WRITE (I W , 1200 ) 

DO  50  I=NL,NH 

50  READ(IR,*)  N,(XIN(I,J),J=1,5) 

NPTS=NPTS+NUM 

IF ( METRO. EQ. 1)  GO  TO  10 

DO  60  I=NL,NH 

XIN (I , 1) =XIN(I , 1) *6894. 757 
60  XIN(I  ,3)=5.*(XIN(I , 3 ) —  32. )/9. 

GO  TO  10 

100  WRITE(IW,  2000) 

READ  (I  R,  * )  NUM,  (IDE(I  )  ,I=1,NUM) 

DO  200  1=1 ,NUM 
II=IDE(I ) 

DO  200  J=1 , 5 
200  XIN  (I  I ,  J)  =0 . 

J=0 

DO  300  1=1 ,NPTS 

IF(XIN(I  ,1)  .  EQ .  0 .  )  GO  TO  300 

J=J+1 

DO  250  K=1 , 5 
XIN ( J  ,  K )  =XIN (I  , K ) 

250  CONTINUE 
300  CONTINUE 

NPTS=NPTS-NUM 
GO  TO  10 

1000  FORMAT ( 2X,48HEDIT  THE  ABOVE  DATA:  ADD,  DELETE,  OR  END? ( A, D,E) ) 
1005  FORMAT ( 1A1 ) 

1100  FORMAT ( 2X, 38HENTER  THE  NUMBER  OF  DATA  POINTS  TO  ADD) 

1200  FORMAT ( 2X , 30HREADY  TO  ADD  DATA  FOR  MATERIAL 
1  /  7X , 21HIN  FORMAT  SHOWN  ABOVE) 

2000  FORMAT ( 2X, 31HENTER  THE  NUMBER  OF  DATA  POINTS 

1  /  2X , 32HTO  DELETE  AND  THE  INDICES  OF  THE 

2  /  11X , 14HROWS  TO  DELETE) 

1150  FORMAT ( 2X, 36HIS  THE  DATA  TO  BE  ADDED  METRIC?(Y ,N) ) 

END 


SUBROUTINE  STORE 

COMMON  DUMMY ( 60 ) , IR, IW, IBUF { 128 ) ,NPTS, XIN(64, 5) ,YD(64) , 

1  COFTO ( 3 ) , TEMPO ,CHISQ ,NPAR, DELTA( 5 ) ,DELMI(5) , 

2  MASK ( 5 ) , GPAR ( 4 ) ,ETPAR(5) ,CNSTR(4) 

DIMENSION  BUF ( 320 ) ,NAME( 3 ) ,BUFA(64) 

DO  100  1=1 ,NPTS 
DO  100  J=l, 5 

100  BUF  (  5*  (I -1 ) +J )  =XIN  (I ,J) 

IBUF(122)=NPTS 
IBUF ( 125 ) =1 
BUFA ( 1 )  =T  EMP0 
BUFA ( 2 )  =C  OFT0 ( 1 ) 

BUFA ( 3 ) =COFT0 ( 2 ) 


i 
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3UFA{ 4  )  =COFTO  (  J ) 

BUFA ( 5 ) =GPAR( 1 ) 

BUFA( 6 ) =GPAR( 2 ) 

BU  FA ( 7 ) =GPAR ( 3 ) 

BUFA ( 8 ) =GPAR( 4 ) 

BUFA ( 9 ) =E  TPAR ( 1 ) 

BUFA ( 10 ) =E  TPAR( 2 ) 

BUFA(11)=ETPAR(3) 

BUFA ( 12 )  =E  TPAR( 4 ) 

BU  FA  ( 1 3  )  =42  TPAR  (  5) 

BUFA ( 14 ) =DELTA ( 1 ) 

BUFA ( 15 ) =DELTA { 2 ) 

BUFA( 16 ) =DELTA ( 3 ) 

BUFA ( 17 ) =DELTA( 4 ) 

BUFA ( 18 ) =DELTA ( 5) 

ICR=15 

LBUF=640 

LIBUF=128 

ICON=3B 

INS=IBUF(121)-1 
IRS=5*I NS 
NAME ( 1 ) =  2HVI 
NAME( 2)=2HBR 
NAME( 3 ) =2H3 

CALL  EXEC (ICR, ICON , IBUF , LIBUF , NAME , INS ) 
NAME( 3 ) =2H4 

CALL  EXEC (ICR, ICON , BUF , LBUF , NAME , IRS ) 
NAME ( 3 ) =2H5 

CALL  EXEC (ICR, ICON , BUFA , LIBUF, NAME , INS ) 
WRITE(IW, 1000 ) 

PAUSE 

RETURN 

1000  FORMAT (2X,22HABOVE  PLOTTED  DATA  AND 
1/2X , 24HPARAMETERS  ARE  STORED  ON 
2/2X, 30HFILES :  VIBR3 ,  VIBR4 ,  AND  VIBR5) 
END 


SUBROUTINE  MODUS 

COMMON  DUMMY(60 ) , IR, IW, IBUF ( 128 ) ,NPTS,  XIN(64, 5) ,YD(64) , 

1  COFT0{ 3 ) , TEMPO ,CHISQ,NPAR,DELTA( 5 ) ,DELMI(5) , 

2  MASK { 5 ) , GPAR ( 4 ) ,ETPAR(5) ,CNSTR(4) 

DIMENSION  A( 5 ) ,ICON(5) ,FR(64)  ,COF(3) 

IY»1HY 

CALL  ERAS 
CALL  HOME 
CALL  ANMOD 
WRITE(I W, 1900 ) 

NFUN*1 

NPAR-4 

TEMP*XIN (1,3) 

FRMIN«ALPHT(T  EMP)*XIN( 1,4) 


r 
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FRMAX=FRMIN 

GMIN=XIN(1,1) 

DO  20  1*1 ,NPTS 
TEMP*XIN(I ,3) 

FRX=ALPHT(TEMP)*XIN(I  ,4) 

IF(FRX.LT. FRMIN)  FRMIN=FRX 
IF  (FRX.GT.FRMAX)  ERMAX=FRX 
IF (XIN(I ,1) .LT.GMIN)  GMIN=XIN(I ,1) 

FR(I  ) =ALOGl (FRX) 

G=XIN(I ,1) 

YD(I )=ALOGl(G) 

20  CONTINUE 

CALL  POLYF(FR,COF) 

A(1)=SQRTIFRMIN*FRMAX) 

Al=ALOGl(  A(l) ) 

A(2)=10.** (COF ( 1 ) +Al* (COF ( 2 ) +Al*COF ( 3 ) ) ) 

A ( 3 )  =.  4 

A(4 ) =.  5*GMIN 

DO  30  1=1 , NPAR 

DELTA (1 ) =0 . 1*  A(I  ) 

IF  (D  ELTA(I  ) .  EQ.  0 .  )  D£LTA(I)=0.1 
30  DELMI  ( I  )  =0 . 01  *D  ELTA  (I  ) 

MASK ( 1 ) =1 
MASK( 2 ) = 1 
MASK( 3 ) =0 
MASK ( 4 ) =1 
WRITE(IW, 2000 ) 

READ ( I R , * )  CNSTR(l) 

WRITE(IW,2100) 

READ(I R, * )  CNSTR( 2 ) 

WRITE (I W, 2200) 

READ(I R, * )  CNSTR( 3 ) 

WRITE dW,  2300 ) 

READ (I R , * )  CNSTR( 4 ) 

CALL  SEARC ( A, NFUN ) 

DO  40  1=1, NPAR 
40  MASK (I ) =0 

50  CALL  SEARC (A, NFUN) 

DO  60  1=1, NPAR 
60  GPAR(I  )  =A (I  ) 

CALL  PLTPK(NFUN) 

WRITE (I W, 1000 ) 

READ(I R, 1005 )  IYN 
IF(IYN.EQ.IY)  RETURN 
DO  70  1=1, NPAR 
70  MASK (I )=0 

WRITE(IW, 1100) 

READ( I R, 1005 )  IYN 
IF(IYN.NE.IY)  GO  TO  300 
WRITE { I W, 1200) 

READdR,*)  NCON ,  (I  CON  (I  )  ,I  =  l,NCON) 

DO  100  1*1, NCON 
100  MASK (I  CON (I ) )=1 
300  WRITE (I W , 1300 ) 

READ(IR,*)  A ( 1 ) 
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WRITE(I W, 1400) 

READ (I R,  * )  A ( 2 ) 

WRITEUW,  1500) 

READ (I R,  * )  A ( 3 ) 

WRITE(IR, 1600) 

READ ( I R ,  * )  A (4) 

GO  TO  50 

1900  FORMAT ( 2X , 16HFIT  MODULUS  DATA) 

2000  FORMAT{ 2X , 41HENTER  THE  MINIMUM  VALUE  OF  FROM  PARAMETER) 

2100  FORMAT ( 2X, 41HENTER  THE  MINIMUM  VALUE  OF  MROM  PARAMETER) 

2200  FORMAT ( 2X, 41HENTER  THE  MAXIMUM  VALUE  OF  MROM  PARAMETER) 

2300  FORMAT ( 2X , 39HENTER  THE  MINIMUM  VALUE  OF  ML  PARAMETER) 

1000  FORMAT( 2X, 25HACCEPT  THESE  VALUES ? (Y , N ) ) 

1005  FORMAT ( 1A1 ) 

1100  FORMAT ( 2X , 34HHOLD  ANY  PARAMETERS  CONSTANT?  (Y  ,N )  ) 

1200  FORMAT ( 2X , 31HENTER  THE  NUMBER  AND  INDICES  OF 
1  /  2X , 23HTHE  CONSTANT  PARAMETERS) 

1300  FORMAT ( 2X , 30HENTER  FROM(Al)  PARAMETER  VALUE) 

1400  FORMAT ( 2X, 30HENTER  MROM (  A2)  PARAMETER  VALUE) 

1500  FORMAT ( 2X ,  27HENTER  N(  A3)  PARAMETER  VALUE) 

1600  FORMAT( 2X , 28HENTER  ML (  A4  )  PARAMETER  VALUE) 

END 


SUBROUTINE  ETA 

COMMON  DUMMY ( 60 ) , IR, IW, IBUF ( 128 ) , NPTS , XIN ( 64 , 5) ,YD(64) , 

1  COFTO ( 3) , TEMPO ,CHISQ,NPAR, DELTA ( 5 ) ,DELMI(5) , 

2  MASK ( 5 ) , GP AR ( 4 ) ,ETPAR(5) ,CNSTR(4) 

DIMENSION  B ( 5 ) ,ICON(5) ,FR(64) ,COF(3) 

IY=1HY 

CALL  ERAS 
CALL  HOME 
CALL  AN MOD 
WRITEdW,  1900) 

NFUN=2 

NPAR=5 

DO  20  1=1, NPTS 
ETAX=XIN ( I ,2) 

YD  (I  )  =ALOGl  (E  TAX) 

TEMP=XIN (I ,3) 

FRX=ALPHT(T  EMP)*XIN(I  ,4) 

20  FR(I  )=ALOG1(fRX) 

CALL  POLYF (FR,COF ) 

B4=-COF(2)/(2.*COF(3) ) 

B ( 1 ) =10 . ** (COF ( 1 ) +B  4* { COF ( 2 ) +B  4*COF ( 3 ) ) ) 

B ( 2 ) = . 35 
B ( 3 ) =- . 45 
B ( 4 ) =10 .  **  B4 
B ( 5 ) =4  . 

DO  30  1=1 ,NPAR 
DELTA  (I  )  =0.1*  fc(I  ) 

I F  (D  ELTA  ( I  )  .  EQ .  0 .  )  DELTA(I)=0.1 
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DELMKI  )»0.01*DELTA(I  ) 

MASK( 1 ) *1 
MASK( 2) =1 
MASK( 3)=1 
MASK ( 4 ) *1 
MASK( 5) =0 
WRITE(IW,2000) 

READ (I R,  * )  CNSTR(l) 

WRITE(IW, 2100 ) 

READ(IR, *)  CNSTR(2) 

WRITE{I W, 2200 ) 

READ ( I R ,  *  )  CNSTR( 3 ) 

WRITE(I W, 2300 ) 

READER,*)  CNSTR(4) 

CALL  SEARC(  B,NFUN) 

DO  40  I«1,NPAR 
40  MASK (I )*0 

50  CALL  SEARC (  8,NFUN) 

DO  60  I»1,NPAR 
60  ETPAR ( I ) *8 ( I ) 

CALL  PLTPK(NFUN) 

WRITE  (I  W, 1000 ) 

READ(IR, 1005)  IYN 
IF(IYN.EQ.IY)  RETURN 
DO  70  I=1,NPAR 
70  MASK (I )=0 

WRITE (I W, 1100 ) 

READ (IRflOOS)  IYN 
IF(IYN.NE.IY)  GO  TO  300 
WRITE (I W, 1200 ) 

READ ( I R , * )  NCON, (ICON(I ) ,1=1, NCON) 

DO  100  1=1, NCON 
100  MASK (I  CON (I ) )=1 
300  WRITE(I W, 1300 ) 

READ(I R, * )  B(  1 ) 

WRITE(I W, 1400) 

READ(IR,  *)  B(2) 

WRITE ( I W, 1500 ) 

READ(I  R,  * )  B(  3 ) 

WRITE (I W, 1600 ) 

READ(IR,  *)  B(  4  ) 

WRITE(I W. 1700 ) 

READ(I  R,  * )  B(  5) 

GO  TO  50 

1900  FORMAT ( 2X, 20HFIT  LOSS  FACTOR  DATA) 

2000  FORMAT { 2X , 44HENTER  THE  MINIMUM  VALUE  OF  ETAFROL  PARAMETER) 

2100  FORMAT(2X,39HENTER  THE  MINIMUM  VALUE  OF  SL  PARAMETER) 

2200  FORMAT(2X,39HENTER  THE  MAXIMUM  VALUE  OF  SH  PARAMETER) 

2300  FORMAT ( 2X , 41HENTER  THE  MINIMUM  VALUE  OF  FROL  PARAMETER) 

1000  FORMAT ( 2X , 25HACCEPT  THESE  VALUES?(Y ,N) ) 

1005  FORMAT (1A1) 

1100  FORMAT (2X,34HHOLD  ANY  PARAMETERS  CONST ANT? ( Y,N ) ) 

1200  FORMAT (2X,31HENTER  THE  NUMBER  AND  INDICES  OF 
1  /  2X, 23HTHE  CONSTANT  PARAMETERS) 

1300  FORMAT(2X,33HENTER  ETAFROL ('Ll)  PARAMETER  VALUE) 
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1400  FORMAT ( 2X , 28HENTER  SL(B2)  PARAMETER  VALUE) 
1500  FORMAT ( 2X , 28 H ENTER  SH(B3)  PARAMETER  VALUE) 
1600  FORMAT { 2X , 30HENTER  FROL(b4)  PARAMETER  VALUE) 
1700  FORMAT( 2X, 27HENTER  C(B5)  PARAMETER  VALUE) 

END 


SUBROUTINE  FCN { X , NFUN ) 

COMMON  DUMMY ( 60 ) , IR, IW, IBUF( 128 ) ,NPTS, XIN( 64 , 5) , YD( 64) , 

1  COFTQ( 3 ) , TEMPO ,CHISQ,NPAR, DELTA (5) ,DELMI ( 5) , 

2  MASK(5) , GPAR(4) ,ETPAR(5) ,CNSTR(4) 

DIMENSION  X( 5) 

F»0. 

DO  30  1*1 ,NPTS 
TEMP=XIN (I ,3) 

FR=ALPHT{TEMP)*XIN(I  ,4) 

DF*YD(I ) -FCT( FR, X,NFUN ) 

30  F«F+DF1t)  F 
CH1SQ=F 
RETURN 
END 


FUNCTION  FCT(FR,X,NFUN) 
DIMENSION  X( 5) 

QQ  TO  (10,20) ,NFUN 
10  FCT»F1(FR,X) 

RETURN 

20  FCT»F2(FR,X) 

RETURN 

END 


FUNCTION  ALPHT  (T  EMP ) 

COMMON  DUMMY ( 60) , IR, IW, IBUF( 128 ) ,NPTS, XIN( 64 , 5) ,YD(64) , 

1  COFTO( 3) , TEMPO ,CHISQ,NPAR, DELTA (5) ,DELMI(5) , 

2  MASK( 5) ,OPAR( 4 ) ,ETPAR(5) ,CNSTR(4) 

ALPHT*10 . *  * ( -1 2 . *  (T  EMP-TEMP0 ) / ( 5  2 5 . /1 . 8+T EMP-T  EMP0 ) ) 

RETURN 

END 
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FUNCTION  F1(FR,A) 

COMMON  DUMMY ( 60 ) , IR, IW, IBUF ( 128 ) ,NPTS, XIN( 64 , 5) ,YD(64) , 

1  COFTO ( 3 ) , TEMPO , CHISQ,NPAR,DELTA(5) ,DELMI(5) , 

2  MASK ( 5 ) ,GPAR{ 4 ) ,ETPAR(5) ,CNSTR(4) 

DIMENSION  A ( 5 ) 

F1=0. 

IF(NPAR.EQ.O)  GO  TO  100 
I F (  A{  1 ) . LT . CNSTR ( 1 ) )  F1  =  1.E10 
IF (  A(  2 ) . LT . CNSTR( 2 ) )  F1=1.E10 
IF (  A(  2 ) . GT. CNSTR ( 3 ) )  F1=1.E10 
I F (  A(  4 ) . LT . CNSTR ( 4 ) )  F1-1.E10 
IF(Fl.NE. 0. )  RETURN 
100  D=l.  +  (  A(1)/FR)**  A(3) 

G=AL0G1 ( A ( 4 ) ) 

Fl=G+( 2 . *  ALOG1 ( A ( 2 )/  A( 4 ) ) )/D 

RETURN 

END 


FUNCTION  F2(FR,B) 

COMMON  DUMMY ( 60 ) , IR, IW, IBUF ( 128 ) ,NPTS , XIN( 64 , 5) ,YD(64) , 

1  COFTO ( 3 ) , TEMPO , CHISQ,NPAR , DELTA ( 5 ) ,DELMI(5) , 

2  MASK ( 5 ) , GPAR{ 4 } ,ETPAR(5) ,CNSTR(4) 

DIMENSION  B ( 5 ) 

F2=0. 

IF(NPAR.EQ.O)  GO  TO  100 
IF(B(1) . LT. CNSTR ( 1 ) )  F2=l.El0 
IF(B{2) . LT . CNSTR ( 2) )  F2=l.ElO 
IF(B (3) .GT.CNSTRl 3) )  F2»1.E10 

IF (B ( 4 ) . LT. CNSTR ( 4 ) )  F2=1.E10 

IF(F2.NE. 0. )  RETURN 
100  A=AL0G1(FR/B(4)  )/B(5) 

C=ALOGl(B(l) ) 

D=1.-SQRT(  l.+A*A) 

F2=C  +  (  (B(2)+B  (3)  )*A  +  (B(2)-B  (3)  )*D  )*B(5)/2. 

RETL.IN 

END 


FUNCTION  ALOGl(X) 
CONST=ALOG(10. ) 

ALOG 1 » ALOG ( X ) /CON  ST 

RETURN 

END 
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SUBROUTINE  VERTL(NUM, LABEL) 
DIMENSION  LABEL(l) 

DO  100  I-1,NUM 

CALL  AOPST ( 1 , LABEL  (I  ) ) 

CALL  MREL(-13,-21) 

100  CONTINUE 
RETURN 
END 


SUBROUTINE  SEARC (X ,NFUN) 

COMMON  DUMMY( 60} , IR, IW, IBUF(128 ) ,NPTS, XIN( 64, 

1  COFTO ( 3 ) , TEMPO, CHISQ,NV,DELTA(5) , DELMI 

2  MASK(5) ,GPAR{4) ,ETPAR(5) ,CNSTR(4) 
DIMENSION  CHI ( 6 ) ,Z(6,5) ,ZBAR(5) ,ZSTAR(5) ,X(5) 


5)  , 
(5) 


ALPHA=1 . 0 
BETA=0 . 5 


GAMMA=  2 . 0 

HUGE*1.0E37 

SIGNI=1.E5 


NVP=NV+1 

NF=0 

1  DO  10  J=1,NV 
10  Z  (NV+1,J)=X(J) 

CALL  FCN ( X ,NFUN ) 

NF=NF+1 

CHI (NV+1 ) =CHISQ 
DO  30  J=1,NV 
DO  20  K=1 i NV 

20  Z  ( J  ,  K )  =X  ( K ) 

SIGN*1 . 0 

23  IF(MASK(J )  .EQ.l)  GO  TO  21 

Z(J  ,  J)=Z  (J  ,  J)+SIGN*D  ELTA(J ) 

21  XS=X(J ) 

X(J)=Z  (J,J) 

CALL  FCN(X,NFUN) 

NF-NF+1 
X(J )=XS 

IF (CHI (NV+1 ) -CHI SQ)  22,30,30 

22  IF  (SIGN  .EQ.  -2.0)  GO  TO  30 

SIGN  =-2.0 

GO  TO  23 

30  CHI ( J ) *CHISQ 

40  JH«1 

JL*1 

DO  80  J*2,NVP 
IF(CHI(J)-CHI(JH) ) 60, 60, 50 
50  JH*J 

60  IF( CHI( J ) -CHI (JL) )70,80,80 

70  JL-J 

80  CONTINUE 

90  DO  120  J*1 ,NV 

ZBAR( J )*0 . 0 


YD( 64) 
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DO  110  K=1,NVP 
IF(K-JH) 100,110,100 
100  ZBAR(J)=ZBAR(J)+Z (K, J) 

110  CONTINUE 

120  ZBAR{ J ) =Z  BAR(J )/FLOAT(NV) 

DO  130  J  =  1 , N V 

IF(MASK( J )  . EQ . 1 )  GO  TO  130 
X( J )  =(  1 . 0+ALPHA ) *  ZBAR ( J ) -ALPHA*  Z ( J H , J ) 
130  ZSTAR{ J ) =X{ J ) 

CALL  FCN ( X ,NFUN ) 

NF=NF+1 

CHIST=CHISQ 

IF  (CHISQ-CHI(JL) >160,180,180 
160  DO  170  J=1 ,NV 

IF ( MASK ( J )  .EQ.l)  GO  TO  170 
X(  J ) =GAMMA*X ( J )  +  ( 1 . 0-GAMMA ) *  ZBAR(  J  ) 

170  CONTINUE 

CALL  FCN  ( X ,NFUN ) 

NF=NF+1 

IF (CHISQ  -  CHIST)  420,230,230 
180  DO  200  J=1,NVP 

IF (J-JH) 190 ,200,190 
190  IF  (CHISQ-CHI  (J))230, 20  0,200 
200  CONTINUE 

IF (CHISQ-CHI ( JH )) 250, 270, 270 
230  DO  240  J=1,NV 
240  Z  ( JH ,  J )  =Z  STAR ( J  ) 

CHI ( JH ) =CHIST 
GO  TO  440 

250  DO  260  J=1,NV 
260  Z  (JH  ,  J )  =X  (J  ) 

CHI(JH)=CHISQ 
270  DO  280  J=1,NV 

IF(MASK(J )  .EQ.  1)  GO  TO  280 
X(J)=BETA*Z(JH,J)  +  (1. 0-BETA)*  ZBAR(J) 

280  CONTINUE 

CALL  FCN ( X ,NFUN ) 

NF=NF+1 

IF (CHISQ— CHI (JH) )4 20, 350, 350 
350  JS=JL 

DO  390  J«1,NVP 
IF (J-JL) 360 , 390 , 360 
360  DO  370  K=1,NV 

Z  ( J  ,  K )  =0 . 5*  (  Z  ( J  ,  K )  +Z  (JL,K) ) 

370  X(K)»Z(J,K) 

CALL  FCN( X , NFUN ) 

NF*NF+1 

IF  (CHISQ-CHI ( JL) ) 38 0,391, 391 
380  JS*J 
391  CHI ( J ) “CHISQ 

390  CONTINUE 
JL=JS 
GO  TO  440 

420  DO  430  J«1,NV 
430  Z(JH,J)«X(J) 
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CHI (JH ) =CH I SQ 
440  DO  470  J=1,NV 

IF (MASK(J )  .EQ.  1)  GO  TO  470 
ZMAX=-HUGE 
ZMIN=HUGE 
DO  450  K=1 , NVP 
ZMAX=AMAX1 ( ZMAX , Z ( K , J ) ) 

450  ZMIN=AMIN1( ZMIN,Z(K, J) ) 

IF  (S  IGNI* ( ZMAX-ZMIN) -AMIN1 ( ABS( ZMAX)  ,ABS( ZMIN) ) ) 470 ,470, 460 
460  IF( ZMAX-ZMIN-ABS (D  ELMI ( J ) ) ) 470 , 470 , 475 
470  CONTINUE 

GO  TO  480 

475  JH=1 
JL=1 

DO  479  J=2,NVP 

IF (CHI (J)-CHI(JH) >477,477,476 

476  JH=J 

477  IF (CHI (J )-CHI ( JL) ) 478, 479, 479 

478  JL— J 

479  CONTINUE 

IF(S IGNI* (CHI (JH)-CHI (JL) )-CHI(JL) )  480,480,90 

480  IF(CHI(JH)-CHI(JL) )490,500,500 

490  JL=JH 

500  CHISQ=CHI  (JL) 

DO  510  J=1,NV 
510  X( J  )  =Z  ( JL ,  J) 

999  RETURN 

END 

SUBROUTINE  GAXES(DELTO) 

COMMON  DUMMY (60) ,IR, IW, IBUF ( 128 ) ,NPTS, XIN( 64 , 5) ,YD(64) , 

1  COFT0( 3) ,TEMP0,CHISQ,NPAR,DELTA(5) ,DELMI( 5) , 

2  MASK (5) ,GPAR( 4 ) ,ETPAR(5) ,CNSTR(4) 

DIMENSION  11(2) ,A(5) ,LABEL(15) 

CALL  ERAS 
ITEN=2H10 
IP=1H+ 

IX=1HX 

NPAR=0 

TEMP=XIN( 1,3) 

FRMIN=ALPHT  (TEMP)*XIN(1,4) 

FRMAX=FRMIN 
GMIN=XIN (1,1) 

GMAX=GMIN 
ETMIN=XIN (1,2) 

DO  50  1=1 ,NPTS 

IF (XIN (I ,1) .LT.GMIN)  GMIN=XIN(I , 1) 

IF (XIN(I ,1) . GT.GMAX)  GMAX=XIN(I ,1) 

IF(XIN(I ,2) .LT.ETMIN)  ETMIN=XIN(I ,2) 

TEMP=XIN(I ,3) 

FR=ALPHT(T  EMP)*XIN(I , 4) 

IFtFR.LT.FRMIN)  FRMIN=FR 
IF(FR.GT.FRMAX)  FRMAX=FR 
50  CONTINUE 

IEFR*ALOGl (FRMIN) 

FRMIN=IEFR-2 
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I EFR-AL0G1 ( FRMAX ) 

FRMAX*IEFR+2 
IGMIN-ALOG1 (GMIN ) 

IGMIN*IGMIN-1 
IGMAX=ALOGl  (GMAX) 

60  I GMAX* I GMAX+ 1 

IDIFG=IGMAX-IGMIN 
I F ( I DI FG . LT . 4 )  GO  TO  60 
GMIN=IGMIN 
GMAX=IGMAX 
IETMI=ALOGl  (E  TMIN) 

ETMIN=IETMI-1 
ETMAX=IETMI-1+IDIFG 
FMIN*0 . 

FMAX=IDIFG 
DELTA ( 1 ) *FRMIN 
DELTA ( 2 ) -FRMAX 
DELTA ( 3 ) *GMIN 
DELTA  ( 4  )  *GMAX 
DELTA  ( 5 )  *ETMIN 
I NDEX* 1 

70  CONTINUE 

IF(INDEX.EQ.l)  CALL  DWIND( FRMIN , FRMAX, GMIN, GMAX) 

I F ( I NDEX . EQ .2)  CALL  DWIND ( FRMIN , FRMAX, ETMIN , ETMAX) 
CALL  TWIND( 140,960 , 50 , 730 ) 

IF (I NDEX. EQ. 2)  GO  TO  117 
FR=FRMIN 

100  DO  105  1*1,10 
XI  =  I 

XFR=F  R+ALOG1 ( XI ) 

CALL  MOVEA(XFR,GMIN) 

CALL  DRAWA ( XFR , GMAX ) 

105  CONTINUE 
FR=FR+1 . 

IF (FR.LT. FRMAX)  GO  TO  100 
G=GMIN 

110  DO  115  1*1,10 
XI=I 

X3*G+ALOGl ( XI ) 

CALL  MOVEA (FRMIN, XG) 

CALL  DRAWA ( FRMAX, XG) 

115  CONTINUE 
G=G+1 . 

IF(G.LT.GMAX)  GO  TO  110 
117  CALL  CHSIZ ( 3) 

DO  120  1*1 ,NPTS 
TEMP*X IN ( I , 3 ) 

FR*ALPHT(T  EMP)*XIN(I ,4) 

G*XIN(I , INDEX) 

EFR*ALOGl( FR) 

EG*ALOGl (G ) 

CALL  MOVEA (EFR, EG) 

CALL  MREL(-4,-6) 

IF (I NDEX. EQ. 1)  CALL  AOPST(l,IP) 

IF (I NDEX. EQ. 2)  CALL  AOPST(l,IX) 


120  CONTINUE 

IF (I NDEX . EQ . 2 )  GO  TO  135 
DO  130  1=1,4 
130  A(I ) =GPAR(I ) 

GO  TO  137 

135  DO  180  1=1,5 
180  A  ( I ) =E  TPAR { I ) 

137  EFR=FRMIN 
140  FR=10.**EFR 

IF  (I  NDEX .  EQ .  1 )  E  G=F  1  (FR,  A) 

I F  (I  NDEX .  EQ .  2 )  EG=F2(FR,A) 

IF  (EFR.  EQ.  FRMIN )  CALL  MOVEA (E  FR,  EG ) 

CALL  DRAWA  (E  FR ,  EG ) 

EFR=EFR+.  05 

IF  (E  FR. LE. FRMAX)  GO  TO  140 
IF (I NDEX. EQ. 1)  GO  TO  145 
CALL  MABS ( 68 , 50 } 

CALL  DRABS (68, 730) 

145  CONTINUE 

I F ( I NDEX . EQ . 1 )  G=GMIN 
I F ( I NDEX . EQ .2)  G=ETMIN 
150  CALL  MOVEA (FRM IN, G) 

IF(INDEX.EQ.l)  GO  TO  155 
CALL  MREL( -68 , 0) 

CALL  DREL( -5 , 0 ) 

155  CALL  MREL(-46,-5) 

CALL  CHSIZ ( 2 ) 

CALL  AOPST ( 2 , ITEN ) 

IEG=G 
CALL  CODE 
WRITE (J ,4000)  I  EG 
CALL  MREL( 1 , 11 ) 

CALL  CHSIZ (3) 

CALL  AOPST ( 2 , J ) 

G=G+1 . 

IF (I NDEX. EQ. 1 . AND. G. LE. GMAX)  GO  TO  150 
IF(INDEX. EQ. 2. AND. G. LE. ETMAX)  GO  TO  150 
I F (I NDEX . EQ . 2 )  GO  TO  170 
MULT*  (TIRMAX- FRMIN  )/10  .  +1 . 

FR-FRMIN 

160  CALL  MOVEA (FR,GM IN) 

CALL  MREL(-26,-25) 

CALL  CHSIZ (2) 

CALL  AOPST ( 2 , ITEN ) 

IEFR*FR 
CALL  CODE 

WRITE( J , 4000 )  I EFR 
CALL  MREL( 1 , 11 ) 

CALL  CHSIZ (3) 

CALL  AOPST ( 2, J) 

FR-FR+MULT 

IF (FR.LE. FRMAX)  GO  TO  160 
170  INDEX* INDEX+1 

IF (I NDEX. LE. 2)  GO  TO  70 

CALL  DWIND( FRMIN, FRMAX, FMIN,FMAX) 
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CALL  TWIND(140,960,50,730) 
CALL  CHSIZ ( 2 ) 

WIDT0-200. 

I F  (D  ELTO .  LT .  20 .  )  WIDT0=100. 

T=TEMPO-WIDTO 

TMAX*TEMPO+WIDTO 

210  FR=ALPHT(T  ) 

EFR=AL0G1(FR) 

CALL  MOVEA  (E  FR,  FMIN) 

EFR=E FR+FMAX 

CALL  DRAWA  (E  FR ,  FMAX ) 

IF  (E  FR.  GE .  FRMAX)  GO  TO  220 
IT=T 

TEM=(T -TEMPO )/50 . 

ITEM=TEM 

IF(ITEM.NE.TEM)  GO  TO  220 
CALL  CODE 
WRITE(I I , 5000 )  IT 
CALL  MREL(— 39 ,3) 

CALL  AOPST( 4,11) 

220  T=T+DELT0 

IF(T  .LE.TMAX)  GO  TO  210 
F=FMIN 

230  CALL  MOVEA (FRMAX, F) 

CALL  MREL( 2,-5) 

CALL  CHSIZ (2) 

CALL  AOPST ( 2 , ITEN ) 

IEF=F 
CALL  CODE 
WRITE ( J  ,  4000 )  IEF 
CALL  MREL( 1,11) 

CALL  CHSIZ (3) 

CALL  AOPST ( 2, J) 

F=F+1. 

IF (F . LE. FMAX)  GO  TO  230 
CALL  CHSIZ (2) 

CALL  MABS( 400,0) 

LABEL ( 1 ) =2HRE 
LABEL ( 2 ) =2HDU 
LABEL ( 3 ) =2HCE 
LABEL ( 4 ) «2HD 
LABEL ( 5 ) =2HFR 
LABEL( 6 ) =2HEQ 
LABEL( 7 ) *2HUE 
LABEL ( 8) *2HNC 
LABEL ( 9) *2HY 
LABEL ( 10 )=2HFR 
LABEL ( 11 ) =2H  H 
LABEL ( 12)S1HZ 
CALL  AOPST (2 3, LABEL) 

CALL  MABS( 430 , 755} 

LABEL( 1)=2HTE 
LABEL ( 2 ) »2HMP 
LABEL ( 3 ) “2HER 
LABEL ( 4 ) *2HAT 
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LABEL ( 5) "2HUR 
LABEL(6)-2HE 
LABEL(7)-2HT 
LABEL ( 8 ) *2HDE 
LABEL  ( 9)  *>2HG. 

LABEL(10)»2H  C 
CALL  AOPST( 20, LABEL) 

CALL  MABS( 5 , 530 } 

LABEL( 1 ) *1HL 
LABEL ( 2) *1H0 
LABEL( 3 )  *1HS 
LABEL ( 4 ) =1HS 
LABEL ( 5) *1H 
LABEL ( 6 ) =1HF 
LABEL( 7 ) *1HA 
LABEL ( 8 ) slHC 
LABEL( 9) =1HT 
LABEL( 10 ) =lHO 
LABEL ( 11 )*1HR 
LABEL( 12 ) =1H 
LABEL ( 13) *1HE 
LABEL(  14 )  s*lHT 
LABEL ( 1 5 ) ®1 HA 
CALL  VERTL( 15, LABEL) 

CALL  MABS( 75 , 635 ) 

LABEL ( 1 ) =1HS 
LABEL ( 2 ) =1HH 
LABEL( 3 ) =1HE 
LABEL ( 4 ) elHA 
LABEL ( 5 ) =1HR 
LABEL ( 6) “lH 
LABEL ( 7 ) *1HG 
LABEL ( 8) *1H 

IF(IBUF(124 ) . EQ. 0)  CALL  VERTL( 8, LABEL) 

LABEL  ( 1 )  *=1HY 

LABEL ( 2)=1H0 

LABEL( 3) *1HU 

LABEL ( 4 ) =1HN 

LABEL ( 5) *1HG 

LABEL( 6 ) *1HS 

LABEL  ( 7 )  =®1H 

LABEL ( 8 ) *1HE 

LABEL ( 9) =1H 

IF (I BUF ( 124 ) .EQ.l)  CALL  VERTL( 9 , LABEL) 

LABEL ( 1 ) *1HM 

LABEL ( 2)®lHO 

LABEL ( 3 ) *1HD 

LABEL ( 4 ) =1HU 

LABEL ( 5 ) =1HL 

LABEL ( 6) slHU 

LABEL ( 7 ) »1HS 

LABEL(8)=1H 

LABEL ( 9 ) “1HN 

LABEL ( 10 ) »1H/ 

LABEL ( 11 ) *1HM 
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LABEL( 12 ) =1H* 

LABEL ( 13 )=1H* 

LABEL ( 14 ) -1H2 
CALL  VERTL( 14, LABEL) 
CALL  MABS( 1008, 520) 
LABEL ( 1 ) =1HF 
LABEL { 2)=1HR 
LABEL( 3 ) =1HE 
LABEL ( 4 ) =1HQ 
LABEL ( 5)=1HU 
LABEL(6)=1HE 
LABEL ( 7 ) =1HN 
LABEL ( 8 ) =1HC 
LABEL( 9) =1HY 
LABEL ( 10) =1H 
LABEL ( 11 )=1HF 
LABEL ( 12 ) =1H 
LABEL ( 13 )=1HH 
LABEL( 14 ) -1HZ 
CALL  VERTL( 14, LABEL) 
CALL  MABS( 0,0) 

CALL  CHSI2 ( 3 ) 

CALL  ANMOD 
PAUSE 
CALL  ERAS 
CALL  RESET 
CALL  CHSIZ (3) 

CALL  HOME 
CALL  ANMOD 
RETURN 

4000  FORMAT  (I  2) 

5000  FORMAT  (I  4) 

END 


SUBROUTINE  PLTPK(NUM) 

COMMON  DUMMY ( 60 ) , IR, IW, IBUF ( 128 ) ,NPTS,  XIN(64,5) ,YD(64) , 

1  COFT0( 3 ) , TEMPO ,CHISQ,NPAR, DELTA( 5) ,DELMI(5) , 

2  MASK ( 5 ) ,GPAR( 4 ) ,ETPAR(5) ,CNSTR(4) 

IF(N  UM-2)  10,20,30 

10  CONTINUE 

WRITE  (I  W,  IQOjO) 

WRITE  (I  W,  1500  )  TEMPO,  (GPAR(I)  ,1=1,4) 

RETURN 

20  CONTINUE 

WRITE  (I  W,  2000  ) 

WRITE  (I  W ,  2500 )  T  EMP0 ,  (ETPAR(I )  ,  1=1 , 5) 

RETURN 

30  CONTINUE 

WRITE(I  W,  3000  )  I  BUF(121)  ,{I  BUFtf  )  ,J=1,30) 

WRITE(I  W,  1000 ) 

WRITE  d  W ,  1500 )  TEMPO,  (G  PAR  (I  )  ,1=1,4) 
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WRITE  (I  W,  2000 ) 

WRITE(I  W,  2500  )  T  EMPO ,  (E TPAR { I }  ,1*1,5) 

WRITE(I  W,  6000 ) 

PAUSE 

RETURN 

1000  FORMAT ( 2X, 47HLOG( M) *LOG( ML)  +  ( 2LOG( MROM/ML) )/ ( 1+(PR0M/FR) * ) 

1  /4X , 4HT0  C, 6X,4HFROM,8X,4HMROM,7X,1HN,9X,2HML, 

2  /15X, 2HA1 , 10X , 2HA2 , 8X , 2HA3 , 8X , 2HA4 ) 

1500  FORMAT ( 2X , F8 . 1, 1PE12. 4, 1PE12. 4, 0PF8. 3,1PE12.4) 

2000  FORMAT ( 2X, 23HAS(L OG(FR) “LOG (FROL ) )/C, 

1  /  2X ,  59H  LOG  <E  TA )  *LOG  (E  TAFROL )  +  ( (S  L4S  H  )A  +  (S  L-S  H )  ( 1-S  QRT  ( 1+A  *  *  2 )  ) )  C/ 
22 , /4X , 4HT0  C , 2X , 7HETAFROL , 4X , 2HSL , 6X , 2HSH , 7X , 4HFROL , 7X , 1HC , 

3  /13X, 2HB1 , 6X , 2HB2, 6X , 2HB3 , 8X , 2HB4 , 8X , 2HB5) 

2500  FORMAT (2X,F8.  1,F8. 3,F8.  3,F8.  3, 1PE12. 4, 0PF8. 3) 

3000  FORMAT(//2X, 15HMATERIAL  CODE  i , 13 ,/2X , 10HMATERIAL  :,30A2) 

6000  FORMAT  ( 2X ,  38H  LOG  <  FR )  =L OG  ( F  )  - 1 2  (T  -T 0 ) /  ( 5  2 5/ 1 . 8+T -T 0 ) // ) 

END 

ENDS 
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